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1.  Introduction 

The  present  study  continues  the  investigation  of  the  author  embodied 
in  NUSC  Technical  Document  TD-7205,  [1].  Here  we  are  again  concerned 
with  underwater  acoustic  scattering  from  the  random  moving  wind-wave 
surface  of  the  ocean.  In  particular,  we  are  concerned  with  high  fre¬ 
quency  0(5-20  kHz),  small  grazing  angle  configurations  0(5-25°),  and  near¬ 
surface  as  well  as  volume  regimes  which  are  bubble-free  vis-a-vis  surface 
effects  [1],  [2]..  Our  efforts  here  are  specifically  directed  to  the 
calculation  of  doppler  shift  and  doppler  spread  produced  by  the  wind- 
driven  ocean  surface.  Included  also  is  a  summary  of  a  more  complete 
model  of  the  soliton  scattering  mechanism  originally  postulated  and 
described  in  [1],  This  more  complete  model  is  based  on  recent  work  of 
Middleton  [3],  [4]  and  Mellen  [4],  It  presents  a  much  fuller  account  of 
the  nonlinear  wind-wave  surface  interactions,  one  of  whose  ultimate 
products  is  the  surface-riding  soliton  ensemble  which  much  empirical 
evidence  suggests  [3],  [4]  provides  the  effective  scattering  mechanism 
results  observed  in  practice  in  "bubble-free"  situations  [2].  A  number 
of  new  calculations  of  backscatter  cross  sections  in  comparison  with 
experimental  data  are  accordingly  included  here,  cf.  [3],  [4]  also,  to 
reinforce  the  plausibility  of  our  postulated  surface  scatter  model,  in 
conjunction  with  the  doppler  analysis,  which  is  the  main  topic  of  this 
report. 


fThis  is  Part  II,  continuing  the  author's  earlier  study,  "Acoustic  Scat¬ 
tering  Cross  Sections  for  Truly  Composite  Wind-Wave  Surfaces:  Scattering 
Without  Bubbles,"  NUSC  Tech.  Doc.  7205,  August  20,  1984  (NUSC,  New 
London,  Conn.  06320). 

*Work  supported  under  contracts  with  the  Naval  Underwater  Systems  Center 
(NUSC),  New  London,  Conn.  06320,  Contracts  N00140-84-M-MS42  and  N00140- 
84-M-NM82. 

**Contractor ,  Physics  and  Applied  Mathematics,  127  E.  91  St.,  New  York, 
NY  10128. 
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1.1  Background 

As  noted  above  (and  in  [1],  [2],  [5],  [6],  for  example),  it  has  been 
widely  observed  for  both  above-surface  (el ectromagnetic-EM)  scattering 
and  underwater  acoustic  scattering  from  wind-wave  ocean  surfaces  that 
such  scattering  is  much  more  intense  0(15-20  db)  in  the  high  frequency 
regime.,  distinguished  by  large  Rayleigh  numbers,  with  moderate  to  fairly 
strong  mean  surface  winds,  U  =  0(5-15  m/sec)  than  classical  theoretical 

a 

scattering  models,  even  those  based  on  composite  surfaces  composed  of 
gravity  and  capillary  components,  can  predict  [1],  [5],  [6].  This  phe¬ 
nomenon  becomes  evident  at  small  grazing  angles,  where  the  short-wavelength 
components  of  the  surface  scatter  dominate.  It  is  especially  noticeable 
in  the  acoustic  cases,  cf.  Sec.  1  of  [1].  For  the  EM  cases,  where  X-band 
radar  is  employed,  the  recent  results  of  de  Loor  [7]  and  Lawner  and  Moore 
[8]  (Fig.  4  especially)  are  similarly  noted. 

In  most  underwater  acoustic  scattering  situations  it  has  been  usually 
postulated  that  a  near-surface  bubble  mechanism,  produced  by  breaking 
waves,  is  responsible  for  the  much  larger  observed  scattering  cross 
sections  at  smal 1  grazing  angles  ([9],  in  addition  to  [5],  [6]).  However, 
recent,  very  precise  experiments  by  Roderick  and  his  group  at  NUSC  [2] 
have  shown  that  these  comparatively  large  acoustic  (backscatter)  cross 
sections  can  occur  without  an  observable  contiguous  subsurface  bubble 
layer,  under  otherwise  comparable  geometrical  and  physical  conditions. 

A  possible  surface  mechanism  which  can  account  for  these  significantly 
larger  scatter  returns  and  doppler  shifts  and  spreads  [2],  when  bubble 
structures  prove  empirically  ignorable,  has  been  suggested  in  [l],  [3],  [4]. 
It  likewise  appears  as  a  critical  product  in  the  generation  of  wind-wave 
surfaces,  as  described  in  [4],  Sec.  2.  The  mechanism  in  question  is  the 
ensemble  of  solitons,  or  hydraulic  bumps  [10],  [ 1 l] ,  which  in  this  model 
are  produced  in  consequence  of  the  nonlinear  wind-wave  interactions  or 
mutual  modulations,  impinging  on  a  thin  wind-driven  water  surface  layer. 

This  latter  is  a  thin  O(centimeter)  moving  surface  layer  which  is  developed 
when  a  nonzero  mean  wind  speed  (U  >0)  is  established,  as  the  fetch  of  the 
surface  becomes  larger;  see  Shemdin  [12];  also  Ramamonjiariosa  [25]. 

[Its  effective  thickness  for  soliton  support  appears  to  be  much  thinner 
o(mi 1 1 imeters ) ,  cf.  Sec.  3  ff.]  Accordingly,  the  scattering  inter¬ 
face  postulated  here  and  in  [1]  consists  of  a  single  wave  surface 
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(usually  artificially  divided  into  so-called  "gravity"  and  "capillary" 
wave  components),  and  a  truly  additional  component,  namely  the  soliton 
"ripples."  These  now  ride  upon  the  moving  wind-driven  surface  portion, 
i.e.,  the  drift  layer,  of  the  gravity-capillary  wave  surface,  in  our 
initial  model  here. 

1.2  Various  Physical  Models 

It  has  been  shown  in  [l]  that  the  various  analytical  methods  and  ap¬ 
proximations  used  in  earlier  analyses  for  the  (scalar)  surface  scatter 
problem  are  essentially  valid  here  and  differ  little  among  themselves  in 
quantitative  results.  The  observed  large  discrepancies  between  theory  and 
experiment,  noted  above,  must  therefore  be  attributable  to  the  choice  of 
physical  model  ([1],  Sec.  1).  Moreover,  on  the  empirical  side,  it  is 
observed  at  once  that  the  radar  data  ([7],  [8]  and  refs,  therein)  indicate 
that  the  level  of  the  wave-number  (intensity)  spectrum  of  the  sea  surface 
in  the  regions  (k  =  5+  to  0.5  rad  cm"*)  depends  markedly  on  the  near¬ 
surface  wind  conditions.  This  implies  that  not  only  U  ,  but  var  U  and 

d  d 

higher  moments,  of  the  turbulent  atmosphere,  become  important  for  a 
proper  description  of  the  sea  surface:  observations  made  at  these  wave¬ 
lengths  will  depend  on  the  particular  sea  surface  actually  generated. 

The  interpretation  of  the  above  is  readily  made:  evidently,  these 

empirical  spectra  are  unsaturated,  i.e.,  are  dependent  on  U„ ,  var  U  , 

.  —  a  a 

etc.,  at  the  shorter  wavelengths.  This  effect  is  evident  in  Fig.  4  of 
[8],  in  Fig.  11  of  Mitsuyasu  and  Honda  [13],  where  wave  gauges  rather 
than  radar  were  employed;  see  also  the  remarks  in  the  more  recent  paper 
[14].  A  common  example  of  a  classical  phenomenological  spectrum  which 
is  saturated  at  the  higher  frequencies  (~u)  )  is  that  of  Pierson  and 
Moskowitz  [15]: 

4  4^4 

2  _5  -0.74gWO:  , 

Wr(f.)  =  Cg  w.3e  ;  C  =  8.10  •  10  J  (indep.  of  U  ,  etc.); 

c,  b  a 

U)s  =  2rfs  (>  0);  Ua>  0.  (1.1) 

Clearly,  this  spectrum  is  saturated,  namely,  independent  of  wind  statis¬ 
tics  (Ua,  etc.)  as  ojs  becomes  large,  for  the  ultragravity  frequencies  on  up. 
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The  empirical  data  noted  above  ([8],  [13],  [14]),  however,  imply 
that  C  -*C(Ua,  var  11^,  ...).  In  particular,  the  data  ([7],  [8],  for 
example)  typically  indicate  spectrum  levels  to  be  0(10-20  db)  above  the 
classical  (~  k  j  average  level,  at  small  angles  (<f>  =  20°-*  10°)  by 
extrapolation.  Accordingly,  this  increased  level,  when  employed  in  the 
classical  theory,  could  possibly  account  for  the  larger  observed  back- 
scatter  cross  sections. 

However,  an  increased  level  alone  cannot  explain  the  observed  doppler 
shifts  and  doppler  spreads  {cf.  Figs.  13,  14  of  [2]  and  Figs.  3.1,  4.1 
following  in  the  present  study}.  Nor  can  such  wave  surface  models,  which 
are  essentially  linear,  account  adequately  for  the  wind-wave  interactions 
which  produce  the  wave  surface.  Thus,  we  propose  here,  as  in  our  earlier 
work  [1],  an  explanation  alternative  to  that  of  simply  raising  the  spec¬ 
tral  level  of  the  gravi ty-capi 1 1  ary  wave  surface  at  the  higher  frequencies. 
The  explanation  chosen  is  that  there  is  an  additional  surface  phenomenon, 
namely,  the  soliton  ensemble,  which  rides  upon  the  moving  drift-layer  por¬ 
tion  of  the  usual  gravity-capillary  wave  surface.  This  ensemble  is  large 
enough  to  account  for  the  observed  backscatter  cross  sections--both  acous¬ 
tic  [1]  and  radar--whenever  a  near-surface  bubble  layer  is  ignorable,  as 
in  Roderick's  experiments  [2].  This  latter  state  can  occur  under  many 
near-surface  wind  conditions,  for  example,  with  large  IT  but  small  var  1)  , 

d  d 

corresponding  to  a  steady  near-surface  wind  with  comparati vely  little 
turbulence,  etc. 

The  reasons  for  choosing  the  soliton  alternative  here  are,  concisely: 
(i).  The  basic  wind-wave  interactions  whereby  wind  energy  is  trans¬ 
ferred  to  the  wave  surface  are  nonl inear,  i.e.,  are  air-water 
"modulation"  processes,  taking  place  primarily  in  the  high- 
wave  number  region  of  the  wave  surface  spectrum,  cf.  [3],  [4], 
and  through  turbulence  transferring  energy  downward  to  the 
smaller  wave  numbers. 

(ii).  Solitons  are  the  result  of  nonlinear  forces,  here  hydrodynamical  ly 
balanced  by  dispersive  forces,  and  they  propagate  nondispersively, 
at  constant  speeds.  These  solitons  have  the  required  "shallow 
channel,"  in  the  form  of  the  thin  drift  layer,  which  is  needed  to 
support  their  generation  and  propagation.  [See  refs.  [6a],  [12], 
[13]  of  [4]  for  further  details,  and  in  particular  [3]  and  [32]. 
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(iii).  Observed  doppler  results  (Figs.  13,  14  of  [2]),  also  Fig.  3.1 
here, are  consistent  with  the  1 i near  frequency  shifts  required 
by  such  nondispersive  elements  as  these  surface  solitons.  (See, 
in  particular.  Fig.  1  of  [25]).  Moreover,  the  observed  doppler 
spreads  [2]  are  likewise  consistent  with  the  soliton  and  drift 
layer  mechanisms,  as  we  shall  show  in  Sections  3  and  4  ff.,  cf. 

Fig.  4.1. 

(iv).  Finally,  the  soliton  model  here  can  also  account  for  the  notice¬ 
ably  larger  than  conventionally  expected  backscatter  cross 
sections,  cf.  Sec.  2,  as  we  have  already  noted  above.  Moreover, 
surface  solitons  can  be  generated  with  long,  intermediate, 
and  short  (spatial)  durations,  as  Olson  et  al  .  [20]  and  others 
[11]  have  demonstrated,  so  that  scattering  contributions  may  be 
expected  over  wavenumbers  0(2tt/25  to  2^70.5  cm"'*')  typically,  in 
the  observed  "high-frequency"  ranges. 

In  all  the  above,  we  emphasize  that  at  this  stage  of  our  experimental 
and  theoretical  knowledge  the  concept  of  the  surface  soliton  ensemble  as 
the  dominant  surface  scatter  mechanism  at  high  frequencies  and  small 
grazing  angles  appears  strongly  plausible,  but  has  not  (yet)  been  proved. 

We  regard  it  as  a  working  hypothesis,  consistent  with  a  broad  spectrum 
of  physical  evidence  [2],  [3],  [4],  [25],  [32]. 

1.3  New  Features  of  the  Present  Study 

Accordingly,  the  new  features  of  this  present  investigation  are, 
conci sely: 

(1) .  A  new  evaluation  of  the  underwater  acoustic  backscatter  cross 

section,  now  as  a  function  of  frequency  o(5-20  kHz),  for 

small  grazing  angles  (9°),  using  the  author's  more  complete 
model  of  the  surface-sol i ton  ensemble  (cf.  Secs.  2  and  9,  ff.); 

(2) .  Calculation  of  the  doppler  shift  5fn,  produced  by  the  combined 

effects  of  a  deep  current,  the  superficial  wind-drift  layer  ([12], 
and  [3],  [4]),  and  the  soliton  movements  on  this  layer  (Secs.  3 
and  7  ff . ) ; 

(3) .  A  determination  of  doppler  spread,  AF,  which  is  seen  to  consist 

primarily  of  spread  due  to  that  in  the  surface  drift  velocity  and 
the  motion  of  the  primary  ocean  surface  (gravity-capillary  surface) 
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and  includes  the  amplitude  modulation  effects  of  the  latter  on 
the  small  scale  components  (Secs.  4  and  10  ff.); 

(4) .  A  comparison  of  the  above  ( 1 ) - ( 3 )  with  recent  experimental  data 

[2],  Figs.  2.2,  3.1,  4.1;  and 

(5) .  A  comparison  of  the  present  theoretical  model  with  "classical" 

results  based  solely  on  the  gravity-capillary  wave  surface,  cf. 
Figs.  2.2,  3.1,  4.1,  and  text. 

The  principal  results  are  presented  and  discussed  in  Part  I,  while  the 
analytical  derivations  are  reserved  to  Part  II,  as  before  [1]. 

1.4  Organization  of  the  Report 

The  Report  is  organized  as  follows:  Section  2  summarizes  our  recent 

results  for  the  (underwater)  acoustic  backscatter  cross  section  as 

l  nc 

a  function  of  frequency,  at  small  grazing  angles,  based  on  the  afore¬ 
mentioned  surface  soliton  model,  which  is  briefly  summarized  here  and 
in  Sec.  7  ff.  Comparisons  are  made  with  Roderick  et  al.'s  recent  data  [2] 
for  bubble-free  regimes  and  Golubin's  earlier  experiments  [16],  which 
represent  different  wind  and  surface  environments.  Section  3  con¬ 
siders  the  theory  and  measurements  [2]  of  doppler  shift,  6fQ,  while  Sec¬ 
tion  4  treats  analogous  results  for  the  associated  doppler  spread,  AF. 
Section  5  completes  the  discussion  of  the  joint  theoretical  and  experi¬ 
mental  results,  which  constitutes  one  of  the  two  main  endeavors  of  this 
study.  The  other  is  the  analysis,  provided  in  Secs.  6-10,  Part  II, 
which  is  based  directly  on  the  results  of  Part  II,  [1],  to  which  the 
reader  is  referred  for  technical  details. 

Part  I Resul ts  and^ Discussion 

2_;^_^IjncjojTej^enjt_^a^cj<scajtt^j^Cxio_s^^^cjtXQJ2LS 

In  the  case  of  underwater  acoustic  backscatter,  the  incoherent  back¬ 
scatter  cross  section  consists  of  three  main  components: 

l  nc 

e(°)  =  jd(°)  +  1  +  &(°)  (2 

inc  |  g-inc  c-incjg  S-inc' 

The  first  two  stem  from  the  gravi ty-capill ary  wave  surface  (G  =  g+c)  and 
have  been  shown  in  [1]  for  the  conditions  of  Roderick's  experiments  [2] 


TR  7635 


7 


to  be  negligible  compared  with  at  high  signal  frequencies  (f  > 

5  kHz)  and  small  grazing  angles  ( 4  =  9°).  The  third  term  represents 

the  backscatter  cross  section  associated  with  our  postulated  soliton 

ensemble  moving  on  the  wave  surface  (G).  Here  6^)  is  a  Bragg  scatter 

(0)  o_inc 

contribution  similar  to  and  proportional  now  to  the  wave  number 

intensity  spectrum  W,;(k|0)<.  of  the  surface-riding  soliton  ensemble. 
Specifically,  we  have 


where 

2 

Rq  =  mean  square  reflection  coefficient  (-  1  here); 

”2 

S  =  mean  square  shadowing  function  (-  1  here,  also); 

kQ  =  Zt</Xq-,  =  wavelength  of  the  (pulsed)  cw  (acoustic) 

excitation,  in  water; 

.  =  "tilt-factor,"  given  by  (7.66a)  and  Sec.  3.1Aof  [1]  . 

bo- 1  DC 


(2.2) 


(2.2a) 


The  soliton  wave  number  spectrum  W^kJO),,  is  found  to  be  ([3],  [4]  and 
Sec.  8.4  ff.) 


,  k  =  2kQ  sin0oT,  (2-3) 


with  0  _  the  angle  of  the  axis  of  the  beam  maximum  with  the  vertical,  so 

that  <j>(=  tt/2  -  0Qj)  is  the  grazing  angle.  The  backscatter  geometry  is 

sketched  in  Fig.  2.1.  Here  £j  =  1  for  semi -i sotropi c  (-tt/2  <  (j>  <  tt/2 ) , 

=  2  for  isotropic  (-tt<  <J><  tt)  excitation  [3]  of  the  solitons  on  the  wave 

2  2 

surface.  (The  quantities  o<.,  are  defined  below  in  Eq.  (2.6a).)  From 

physical  considerations  we  usually  reject  the  isotropic  model,  since 

the  soliton  surface  moves  downwind.) 

(0)  ~(0) 

The  tilt-factor  11^'  (=  16  N^;)  accounts  for  the  amplitude  modulation 


2  2 
TT  £  t  Qq  Oi 

w  (k|o)  = - ^ 

i  5  [i  +  {kajzrr 
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Figure  2.1  Backscatter  Geometry 


of  the  small-scale  surface  waves  by  the  gravity  wave  component  of  the 
surface  and  is  for  these  backscatter  geometries: 


“gs*  1  16{3<4s'in',0oT  +  6aLs1"2eoTcos20oT  +  cos\t>’ 


with  4»oT  =  ti/2;  oGx  =  oGy. 


(2.4) 


Here  the  rms  large-scale  slopes  are  found  numerically  from  Cox  and  Flunk's 
result  [17] 

=  (3.0  +  5.12  U  )  • 10"3.  (2.5) 

d 

The  reflection  coefficient  is  essentially  unity  here,  and  the  shadowing 

2 

coefficient  is  likewise,  e.g.,  S  -  1,  as  long  as*  0Qj<  85°  ( cp  >  5  °  )■ . 

*For  cases  where  £  <  1  and  thus  must  be  considered,  see  Bass  and  Fuks 

[18] ,  Fig.  59,  and  recent  work  on  underice  scattering.  Sec.  3.2C  of 

[19] . 


JG(x=y) 


,  (K 

3x 
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Combining  (2.4)  and  (2.3)  in  (2.2)  gives  the  desired  backscatter  cross 
section  explicitly  (with  =  1  now) 


9(0) 


2  2/p 

0  {\  A  p  p  p  44 

-  (3aGxsnn  9oT+6oGxsln  e0jcos  e0j  +  cos  3oT^ko 


TT 


[l+(ahkos1r,0oT)2]2 


(2.6) 


where 


=  mean  square  intensity  of  the  surface  soliton  ensemble  (cf.  Sec.  7.3); 


=  mean  square  "wave"  length  of  the  solitons. 

Results  of  recent  experiments  and  Eq.  (2.6)  of  the  present  theory 
for  acoustic  backscatter  are  presented  in  Fig.  2.2,  for  0Qj  =  81°  (or 
<$>_=  9°_±  and  N^  (=  N^/16)  =  1.68  10”2  for  lT-  =  10  m/sec.  Here 

-  52  =  1  and  +2.5  db  has  been  added  for  the  Kirchoff  correction  at 
these  small  grazing  angles  (all  curves),  cf.  Table  3.1  and  "footnote" 
thereto  in  [1].  The  fit  of  theory  to  experiment  is  made  for  the  above 
geometry  over  a  range  of  5-20  kHz.  Narrow  band  signals,  i.e.,  compara¬ 
tively  long  cw  pulses,  are  used  in  transmission  and  reception  by  Roderick 
et  al . ,  cf.  Fig.  10  [2].  Explosive  charges  were  employed  by  Galubin  [16], 
also  Fig.  1.26  of  [9],  with  narrow-band  reception.  In  Figure  2.2  the  O 
represent  the  Roderick  data  [2],  Fig.  10,  while  the  ®-points  are  Galubin's 
data  [16].  Included  in  Fig.  2.2  are  some  typical  theoretical  backscatter 
cross  sections  based  on  purely  capillary  resonance  returns  ([9],  Sec. 

9.11,  Figs.  1.26,  9.10,  and  [1],  Eq.  (3.18)). 


(2.6a) 


2.1  Discussion  and  Results 

Clearly  evident  are  the  different  near-surface  wind  states  for  the 
two  data  sets.  Thus,  although  the  mean  wind  speeds  U~  were  essentially 

d 

the  same  (  20  knots  =  10  m/sec)  in  both  instances,  different  wind  states 
were  in  force,  as  quantified  by  different  values  of  var  U  (and  higher 

d 

moments)  although  not  recorded  at  the  time.  We  infer  this  from  the  dif¬ 
ferent  thickness  of  the  wind-drift  layers  (-o^)  and  hence  average  length 
of  the  surface  solitons,  and  by  the  different  rms  soliton  amplitudes 
(~as).  As  Shemdin  has  shown  ([12],  Figs.  1  and  2),  the  greater  the  wind 
speed  U  (and  U_),  the  faster  movinq  but  thinner  is  the  resultant  surface 

a  a 
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Figure  2.2  Incoherent  Backscatter  Cross  Sections  (Middleton  Soliton  Model,  1984-,  [3],  [4]) 
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drift  layer,  and  for  fetches  sufficiently  large  the  drift  layer  velocity 
becomes  nearly  constant.  This  thinner  layer  can  support  larger  and  faster 
solitons,  depending,  of  course,  on  the  intensity  and  velocity  of  the 
initiating  atmospheric  shock  waves  [3],  [4],  [32]. 

Accordingly,  on  the  basis  of  our  postulated  soliton  model  here  we 

can  deduce  from  Fig.  2.2  that  the  Roderick  data  are  associated  with  a 

rather  noticeable  secular  spread  of  wind  speeds  (~  large  a^).  To  achieve 

the  observed  backscatter  cross-sectional  values,  an  appropriately  larger 

2 

soliton  intensity  (cr<.)  is  then  required.  With  the  present  model  this 
means  that  large  but  comparatively  uniform  shock  wave  intensities  are 
maintained,  leading  to  "large"  soliton  velocities  (v  )  with  small  vari¬ 
ances,  at  levels  consistent  with  the  data  [2].  On  the  other  hand, 
Galubin's  data*  imply  a  narrower  spread  (o^)  in  the  drift  layer,  which 
suggests  smaller  var  U  .  In  addition,  smaller  values  of  v  are  needed, 
producing  smaller  levels  of  soliton  intensity  (a^)  for  nearly  the  same 
observed  backscatter  cross  section  at  the  higher  frequencies.**  From 
this  we  conclude  that  different  near  surface  environmental  states  are 
represented  by  the  two  sets  of  data,  even  though  U  is  effectively  the 

d 

same.  This  is  consistent  with  our  observation  in  Sec.  1.2,  regarding 
the  unsaturated  Wave  spectra  as  frequency  is  increased. 

In  general,  we  also  note  that: 

(i).  Both  the  theoretical  and  empirical  (incoherent)  backscatter  cross 
sections  become  independent  of  frequency  as  kQ-*<»,  cf.  (2.6). 

(ii).  The  potential  existence  of  "long"  solitons  0(15-25  cm),  as  well 
as  shorter  ones,  supports  the  possibility  of  backscatter  returns 
at  the  lower  "high"  frequencies  0(5  kHz),  cf.  Fig.  2.2,  as  well 
as  at  20  kHz,  etc.  For  examples  of  such  solitons  produced  in 
tank  experiments,  see  Figs.  2-6  of  Olson  et  al .  [20],  esp.  Fig. 

2  where  =  5.0  cm  soliton  length,  Lq  -  18  cm,  and  height  ~ 

3  cm,  is  typical  of  solitons  capable  of  Bragg  scatter  at  ~5  kHz. 


*These  are  regarded  as  bubble-free  here:  Galubin's  experimental  procedure 
indicated  negligible  bubble  effects  in  the  volume  but  was  not  capable  of 
discriminating  a  possible  bubble  surface  layer. 

**The  noticeably  larger  value  of  at  5  kHz  may  be  attributable  to  the 
growing  contribution  of  the  gravity-wave  "facet"  terms  for  these  rather 
small  Rayleigh  numbers. 
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2.2  Remarks  on  the  "Bubble-Free"  State 

Finally,  a  number  of  observations  regarding  the  "bubble-fre"  nature 
of  Roderick's  experiments  [2]  are  in  order: 

(1) .  The  resolution  distance  is  where  xp  =  0.4  msec  is  the 

pulse  length  used,  directed  normally  to  the  wave  surface  from 

below,  so  that  x  c  /2  =  0.30  m  is  the  resolution  distance  here 
P  o 

(cf.  Fig.  3  of  [2]  and  corrected  text).  Outside  this  distance 
about  the  mean  surface  elevation  (c  =  0)  vis-a-vis  the  projector 
no  bubbles  distinguishable  from  the  much  lower  volume  reverbera¬ 
tion  were  observed. 

More  precisely,  operating  at  vertical  incidence  (0  y  =  0) 
the  receiving  system  was  self-noise-limited  such  that  if  bubbles 
were  present  in  the  volume  below  the  water/air  surface,  their 
backscatteri ng  strength  would  have  to  be  more  than  -47  db  to  be 
apparent.  This  compares  with  -30  db  actually  observed  for 
surface  scatter  (>  15  kHz),  cf.  Fig.  2.2.  No  volume  scatter  at 
or  above  -47  db  was  noted  as  a  precursor  to  the  surface  return, 
which  indicates  that  the  empirical  backscatter  returns  in  this 
experiment  were  well  dominated  by  surface  scattering  mechanisms. 

(2) .  Moreover,  if  near-surface  bubbles  are  present  in  quantities  suf¬ 

ficient  to  mask  the  surface  and  provide  sufficient  return  to 
account  for  the  observed  levels,  at  these  high  frequencies  and 
small  grazing  angles,  the  work  of  Glotov  [21]  and  Thorpe  [22] 
has  shown  that  these  bubbles  will  appear  in  variable  layers  0(2-4 
meters)  thick  immediately  below  the  surface.  In  particular,  see 
Fig.  5D,  p.  165  of  [22].  This  is  much  deeper  than  the  0.3  m 
resolution  distance  achieved  by  Roderick  et  al.,  cf.  (1)  above: 
no  such  layers  were  observed  [2]. 

(3) .  When  a  masking  bubble  layer  is  present,  there  will  be  little  or 

no  observable  doppler  shift,  neither  for  the  now-masked  resonant 
surface  wave  components  nor  any  surface  drift  currents.  The  data  of 
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[2],  cf.  Figs.  11,  13,  show  a  very  distinct  and  measureable  dop- 
pler  shift,  in  addition  to  an  expected  doppler  broadening  of  the 
spectrum,  Fig.  14,  [2]. 

Thus,  as  we  have  noted  above,  and  in  [1],  [3],  [4],  [32]  the  conclusion  is 
drawn  that  some  specific  surface  mechanism,  here  the  proposed  soli  ton 
ensemble  (from  the  additional  argument  presented  in  Sec.  1.2),  is  re¬ 
sponsible  for  the  observed  scatter  results,  when  bubbles  are  ignorable. 

We  emphasize  that  the  surface  soliton  model  advanced  here,  and  in  [l]-[4], 
[32],  is  so  far  inferential ly  supported  by  the  data  and  physical  argument, 
but  has  not  yet  been  directly  proved.  A  direct,  quantitative  demonstra¬ 
tion  of  the  soliton  wave-surface  component  remains  to  be  demonstrated 
explicitly,  in  the  process  of  scattering.  Nevertheless,  as  we  show 
here,  our  soliton  surface  model  is  able  to  account  for,  and  predict,  not 
only  cross-sectional  values,  but  doppler  shifts  and  spreads  where  the 
classical  models  fail,  vide  Table  5.1  ff. 

3 .  Xhe,  Me_a  n_  D  0  o  ^j^^JTheorv^^an  d__E  x_p  ej^me  njt 

From  the  analysis  of  Section  7,  for  the  geometry  of  Fig.  7.1  which 
applies  for  the  experiments  of  Roderick  et  al .  [2],  we  have  the  (down) 
doppler  shift  (7.18a): 


2sin6 


o 


A, 


(3.1) 


o 


whose  structure  and  derivation  are  discussed  (7.18)  ff. 
The  measured  parameters  of  the  model  are  [2]: 


vc  =  (mean)  deep  current 
<j>  =  angle  of  current 

4,  =  mean  wind  direction 


l)  =  mean  (near  surface)  wind  speed  =  10  m/sec. 

a 


W 


=  26.0  cm/sec 


51.5°;  sin4>c  =  0.783; 

116°  ;  si n$  =  0.899; 
w 


(3.2a) 
(3.2b) 
(3.2c) 
(3. 2d) 


In  addition,  the  related  backscatter  geometry  has 


<i>oT  =  7T/2  ;  eoT  =  81  (<}>  =  9°  grazing  angle). 


(3.3) 


TR  7635 


14 


Furthermore,  we  assume  that  the  surface  drift  layer  moves  in  the 
mean  wind  direction  ($),  and  that  the  measurements  at  f  =  20  kHz  are 
probably  the  most  accurate.  From  the  ensemble  of  measured  <5f  's,  repre¬ 
sented  by  the  points  ®  in  Fig.  3.1  (taken  from  Fig.  13  of  [2]),  we  find 
that 


<{fo>expt.  ■  12-4  Hz’  <V7-5cn,), 


(3.4) 


and  since 

l<5f 


'2sin0 


oT' 


X 


oD’ 


oD 


xpt 


47.0  cm/sec 

=  cc  +  v  s i n 4>  +  v  .sin4> 

S  c  c  d  w 


(3.5) 


where  sin  81°  =  0.988. 

Our  next  step  requires  an  assessment  of  v^.  For  this  we  consider 
a  range  of  reasonable  values,  in  part  suggested  by  the  "rule-of-thumb" 
that  the  mean  drift  speed  is  0(1-4%)  of  the  mean  wind  speed  l)  (Bass  et 

_  a 

al .  [29];  [25],  Fig.  5).  Having  specified  a  value  of  v^  we  can  then 

use  (3.5)  to  obtain  an  estimate  of  the  (total)  soliton  speed,  c<-  = 

c  +  v  ,  cf.  Eq.  9.2  et  seq.  From  the  backscatter  cross-sectional  data  of 
oo  o  v 

Fig.  2.2  we  find  that  the  rms  soliton  wavelength  =  (|_  ) 2  =  5.0  cm. 

From  this  and  the  estimate  of  c<~  we  then  can  show  that  the  incremental 

soliton  velocity  v  ,  due  to  the  nonlinear  mechanism  involved  [cf.  Sec.  9] 

is  ignorable  vis-a-vis  c^  =  c  ,  and  so  deduce  directly  the  effective 

thickness  of  the  drift  layer,  h  supporting  the  soliton  mechanism. 

We  begin  with  the  deep  current  and  obtain 

v^sin  51.5°  =  26  sin  51.5°  =  26(0.783)  =  20.3  cm/sec.  (3.6) 

Table  3.1  shows  values  of  c<~,  hg^^  and  vQ,  for  a  range  of  reasonable 
values  of  v^:  the  shaded  portions  of  Table  3.1  appear  to  be  appropri¬ 
ate  to  the  model  and  circumstances  described  here.  In  addition,  the 
relations  used  to  determine  c<-,  ^-f,  v0’  are  (from  Sec.  9  ff.) 

CS  =  S  +  v0  =  /9heff  +  V  vo  =  4  heff  /9/3Lo  ;  9  =  980  cm/sec2’ 

(3.7) 
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igure  3.1  Doppler  shift  -<5f  ,  as  a  function  of  frequency,  for  the  soliton 

(Eq.  13.1)  and  classical  models  (3.8),  including  surface  drift  and 


deep  current  speeds.  Theory  and  experiment  [2]. 
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Table  3.1: 


(U  =  10  m/sec) 

a 


vd 

v j=v  ,sin$ 
d  d  w 

CS  voD'20-3 vd 

hefflc(/9 

n6'4 

eff 

vo 

10  cm/sec 

9.0  cin/sec 

17.7  ctn/sec 

3.2  ran 

5.79- 10"2 

-2 

9.7*10  cm/sec 

15.0 

iil3.5 

K 

13.2 

1.8 

1.37- 10"2 

2.3- 10'2 

20.0 

:  18.0 

8.7 

0.77 

1.6  -10"3 

2.7* lO-3 

£5.0 

j 

22.5 

4.2 

0.018 

4.3  -10"5 

-5 

i  7.2-10  3 

28.0 

25.2 

1.5 

0.0023 

-7 

2.5  -10 

4.2-10"7 

and  vQ  «  cq  (-  c,.)  here,  with  LQ  =  21.^  =  ch  =  5.0  cm  [cf.  Fig.  2.2  here]  from 
the  fit  of  the  backscatter  cross-sectional  data  of  Roderick  et  al .  [2].  Rea¬ 
sonable  surface  layer  thicknesses  and  speeds  of  the  surface  soliton  ensemble 
are  noted  in  the  shaded  regions  of  Table  3.1,  as  well.  [See,  however,  footnote, 
p.  50.] 

If  instead  of  the  soliton  surface  model  we  employ  the  purely  classi¬ 
cal,  linear  model  of  a  gravity-capillary  wave  surface,  we  find  that  in 
place  of  (3.1)  the  doppler  shift  6 f Q  is  now  given  by  (7.22b),  namely, 


Sf 


ol grav 


Ig  sin9  ,  2sin6  ,  _  _  i 

^r; —  +  ^r —  tv1n*c  *  VinV: 


(3.8) 


where  >  5  cm,  so  that  gravity  waves  predominate.  Equation  (3.8)  in¬ 
cludes  the  contributions  of  the  observed  deep  current,  and  the  surface 


drift.  Figure  3.1  shows  6f  for  various  mean  drift  speeds 

o  y  rd v 

0<  v^  <20  cm/sec . 

From  the  results  above,  namely  (3.1)  with  (3.5)  and  Table  3.1,  and 
(3.8),  we  see  that  the  soliton  surface  model,  which  implies  a  1  i near 
doppler  shift,  fits  the  data  well  above  10  kHz  and  is  not  excessively 
off  at  the  lower  frequencies.*  On  the  other  hand,  even  for  zero  drift 


*As  before,  cf.  Fig.  2.2  (Galubin's  data),  the  behavior  at  5  kHz  departs 
from  theory  more  noticeably  than  do  the  data  at  other  (higher)  frequencies 
We  suspect  that  this  is  the  result  of  the  small  Rayleigh  numbers  involved, 
which  reflect  the  fact  that  now  contributions  from  the  large-scale  surface 
components  are  not  entirely  ignorable. 
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speeds  the  "classical"  cases  do  not  agree  nearly  as  well  with  the  data, 
and  for  larger  drift  speeds  the  discrepancies  are  quite  pronounced.  We 
can  say,  at  least,  that  the  soliton  surface  model,  which  is  nondi spersi ve , 
(cf.  Sec.  1,  Sec.  8.3)  and  hence  exhibits  disturbances  moving  at  a  con¬ 
stant  speed,  is  consistent  with  the  careful  observations*  of  Roderick 
et  al .  [2],  with  respect  to  (mean)  doppler  shift.  This  is  in  addition 
to  the  plausibility  of  the  model  demonstrated  by  the  backscatter  data, 
as  discussed  above  in  Section  2. 


4.  Doppler  Spread,  AF:  Theory  and  Experiment 

The  doppler  spread,  measured  between  10  log^e-*  points  below  the 
spectral  maximum  in  accordance  with  Eq.  (10.10)  and  following  the  em¬ 
pirical  procedure  of  Roderick  et  al .  ([2],  Fig.  14),  is  found  from  the 
analysis  of  Sec.  10  ff.  (Part  II)  to  be 


log  C(=10  logine) 

A  F  =  - — - ( 2a  k  a . )  ^  1  + 

tt  v  oz  o  t,  1 


which  can  be  more  compactly  written 


3. 


R, 


(2“oza;/Co)2  ko(2aozaj)2 


AF  =  Ak  1  +  B  +  C/k~" 
o  *  o 


(=  \|A2k2  +  A2Bk2  +  A2C  ), 
u  0  0 


(4.1) 


(4.2) 


where 


n  _  V2  log  C  ,n  ^  1.714  x  , 

A  =  - — 3 —  (2a  a;)  =  (2a  a*  k 

it  oz  4 7  tt  v  oz  z, '  1 


a-  =  Eq.  (8.14);  aQZ  =  cos0qT 


2  ill  \  2 

fg/R-  =  ^  11  02  c  '  ' 


B  =  aQ/R.  =  //  — 1  ;  where 


03  =  °c  +  8d’  Eq-  (8>4); 


C  =  RG/(2aozar)  ;  RG  =  (Eq-  8‘29b)  ■ 


(4.2a) 

(4.2b) 

(4.2c) 


Here  the  term  (~AkQ)  represents  the  phase  modulation  (of  the  carrier 
frequency  f  )  which  is  produced  by  the  random  movement  of  the  large-scale 
(i.e.,  gravity-capillary)  wave  surface,  on  which  the  soliton  ensemble 

p 

and  drift  layer  move.  The  second  term,  involving  A  B,  embodies  the  phase 

p 

modulation  generated  by  the  random  variations  (~a£)  in  the  deep  current 


*See  footnote,  page  16. 
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2  2 

speed  and  in  the  drift-current  speed  (~cr  ).  Finally,  the  term  A  C  gives 
the  contribution  of  the  amplitude  modulation  created  by  the  random  tilt¬ 
ing  of  the  local  large-scale  surface  on  which  the  small-scale  components 
ride  and  move.  Unlike  the  phase  modulations  here,  this  amplitude  modu¬ 
lation  is  independent  of  the  incident  wave's  frequency,  cf.  (4.2). 

2 

To  determine  and  Rg  two  spectral  spread  measurements  are  required. 
(In  principle,  we  could  calculate  Rg  if  we  knew  how  to  truncate  the  point 
spectrum  of  the  surface  elevation  Cr,  cf.  (8.27a),  but  this  does  not 
appear  generally  possible.)  Thus,  both  and  Rg  represent  truly  empiri¬ 
cal,  albeit  physical  quantities,  to  be  found  in  practice  by  measurement, 
as  we  shall  do  here.  Accordingly,  if  we  select  two  different  measure¬ 
ments  of  doppler  spread  AF^,  AF^,  for  the  given  conditions  of  the  experi¬ 
ment:  U  .  a  ,  k  ,  (a*)  appropriate  to  "high"  frequencies  f  >0(5  kHz) 
and  small  grazing  angles  (©^>70°),  we  can  write  from  (4.2) 


Vkol  VI  +  B  +  C/&  ;  A2  =  k02  lf1  +  BtC/k02  ;  Al,2 


Solving  for  B  and  C  is  routine.  The  results  are 


AFij2/A. 

(4-3) 


B  - 


.2  1 2 
ko2  ‘  kol 


1  (>  0); 


aV  -aV 

c  -  *  f  p1  (>0). 

-  v 


(4.4) 


since  from  (4.2b,c)  B,  C  must  be  positive.  We  give  several  examples  of 
the  explicit  calculations  in  Sec.  4.1  ff. 


4.1  General  Model  I,  Eqs.  (4.1),  (4.2):  (Roderick's  Experiments  [2]) 

For  the  specific  parameters  of  the  recent  experiments  of  Roderick 
et  al .  [2]  we  have 

©of  =  81°  (.*•  4)  =  9°  grazing  angle);  mean  wind  speed  Ua  =  10  m/sec;' 

cQ  =  1.5  •  103  m/sec. 

a-  =  2.09  -D72  •  10-3  =  0.209  (PM  Spectrum,  cf.  Eq.  (8.14)); 

Q  a 

a*  =  0.457; 


a 


oz 


=  cos0oT 


0.157-  ; 


a2  =  5.42  •  10_2,  Eq.  (8.10);  /.a  =  0.233; 

X  X 

(a|)G  =  1.42  *  10"5  (T4  =  0.142  (m2);  /.  (a 


0.377. 


(4.5a) 

(4.5b) 
(4.5c) 
(4 . 5d ) 
(4. 5e) 
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•’•A  =  (0.313)  (0.457)  = 


7.82  -10 


B  =  a2/ (9.54  •  10"8) 


C  =  Rr/ (2.05  •  10"2) 


(4.6) 


The  elements  of  the  calculations  based  on  (4.4)  are  given  in  Table  4.1: 


Table  4.1  Calculated^ ^op^l^r^Spreads AF (  ^ 

I. 


^ - Expt'l  >Curve(s)  Fig,  4.1 


f 

0 

AF  (Xpt'l ;[2],  Fig.  14) 

k 

0 

7T 

j  O  IX) 

1 

A  = 

<V  A 

A(f  ) 

0 

■’'AF10,20 

AF5,20 

20  kHz 

23.7  Hz 

83.78 

7019 

3.03* 102 

9.185- 104 

A 

23.7  Hz 

23.7  Hz 

15 

22.0 

62.84 

3949 

2.81- 10* 

7.915*  nr 

(20.6) 

(20.4) 

10 

18.0 

41.89 

1755 

2.30- 102 

5.298* 104 

18.0 

(17.4) 

5 

15.0 

20.95 

439 

1.92- 102 

3.679- 104 

(16.3)  . 

15.0 

o-  calculated  from  (4.3),  based  on  AF^,  AF2o’  or  AF5S  AF20’  aS  indicatecl* 


Table  4.1a  Auxiliary  Parameters,  Eq.  (4.4) 


B10,20 

10C20 

B5,20 

5C20 

6.38 

4.00  104 

8.37 

3.32  104 

*'•  ae 
••  rg 

7.83* 10~4rad 

8.20* 102sec-2 

8.94* 10~4rad 

6.81- 102sec"2 

AFah  =  Eq.  (4.8) 

AFe  5  f„°B 

15.7  Hz 

15.6  Hz 

17.9  Hz 

14.2  Hz 

The  results  of  applying  (4.3)  are  shown  in  Fig.  4.1;  the  solid  curve 
(for  I  =  Model  I)  is  based  on  the  selection  of  the  two  data  points  AF^g  = 
23.7  Hz  and  AF^g  =  18.0  Hz,  where  the  calculated  points,  following  this 
choice  of  AF^g,  AF^g,  are  indicated  by  the  circled  numbers  in  Table  4.1. 
Similarly,  the  calculated  points  in  the  dotted  curve (Model  I)  are  obtained 
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Figure  4.1  Doppler  spread  AF,  theory  and  experiment,  Eq.  (4.1),-  (4.2), 
for  the  soliton  model.  Classical  (gravity-wave)  and  other 
results  (4.9)  and  Table  4.2. 
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from  the  selection  of  AFon,  AFC.  The  fit  appears  better  for  the  first 
choice  (solid  line);  our  selection  of  AF^q,  in  any  case,  is  reasonable, 
as  this  result  is  probably  more  accurate  than  the  others,  both  from  the 
experiment  and  from  the  fact  that  f  is  larger  and  therefore  the  ob¬ 
served  AF^q)  is  more  accurate. 

The  auxiliary  quantities  B,  C  are  shown  in  Table  4.1a,  along  with 
estimated  values  of  and  R^.  If  we  define  a  measure  of  spectral  spread 
(not  the  same  as  used  in  (4.1)  et  seq.,  cf.  Eq.  (10.10),  but  proportional 
to  it)  by 


AF3  =  f0ae,  cf.  (8.3)  ,  (4.7) 

we  see  (Table  4.1a)  that  at  f  =  20  kHz  the  rms  spread  due  to  the  com¬ 
bined  effects  of  deep  current  and  dr.ift  velocity  variation  is  o(16  Hz). 
Similarly,  we  see  that  if  the  spread  due  to  the  AM  effects  of  the  random 
tilt  is  defined  by 


AF 


AM 


_  1  im 


AF  =  A  /C  = 


1.714 

TT 


/R, 


(4.8) 


then  the  resulting  doppler  spread  due  to  this  mechanism  is  o(14,  15  Hz), 
cf.  Table  4.1a.  [The  overall  resultant  doppler  spread,  as  given  by 
(4.1),  (4.1),  of  course,  is  not  simply  the  sum  of  the  individual  spreads, 
because  of  the  way  the  various  spreading  mechanisms  interact  physically, 
vide  Secs .  8-10.  ] 

As  frequency  becomes  less,  our  results  (4.1)  approach  a  nonvanishing 

value.  This,  of  course,  does  not  actually  occur  physically  below  some 

f  ~0(+2  kHz),  where  the  dominant  Bragg  scatter  term  is  now  that  of  the 

large-scale  component  (not  evaluated  here).  Thus,  our  results  for  Model  I 

here  are  not  to  be  extrapolated  to  frequencies  below  o(l-2  kHz).  On  the 

other  hand,  for  frequencies  much  above  20  kHz  the  spatial  terms  in  F0  „  , 

2 -a 

(8.8),  appear  as  they  do  in  the  tilt  factor,  K  /g)(Ar,x),  cf.  (7.5a).  We  no 
longer  have  a  simple  soliton  "wave"  surface,  but  one  (cf.  Z  in  (7.8)  et 
seq.)  where  the  dispersion  law  K  is  no  longer  simple.  This  changes  the 
doppler  shift,  6f  ,  cf.  Sec.  7.  (vs.  Secs.  7.2,  7.3),  but  does  not 

u  In)  .  - 

change  the  x-structure  of  mrT  (Tlf0)>  cf.  ( 10 . 3 )  -  ( 10 . 7 ) ,  where  Ks+  K-, 
now:  the  scale  °f  fvx-<x>  is  altered’  but  not  the  spreading  terms 
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2 

(~  o^),  (10.2).  We  conclude,  therefore,  that  as  frequency  (f  )  is  in¬ 
creased,  i.e.,  as  k  -*-00,  Equation  (4.1)  (and  (4.2))  still  provide  the 

0  -j  • 

measure  of  doppler  spread,  i.e.,  |<Q^00  AF-*AkQ  /l+B  ~  0.212  kQ,  cf. 

Fig.  4.1.  Moreover,  at  very  high  frequencies,  the  backscatter  cross 

section  ,  (2 .6),  becomes  independent  of  frequency,  [whereas,  as  k  -*0+, 

JO)  4  5  0 

°S  ~  ^o’  a"*^elt  that  this  effect  is  overwhelmed  by  an  analogous  Bragg 

scatter  from  the  gravity-capillary  component.] 

Finally,  Fig.  4.1  also  shows  our  Models  II,  III,  where  we  have 

respectively  dropped  the  contributions  (~  B,  C)  and  (~  C),  in  (4.2). 

Table  4.2  lists  pertinent  values  of  AF  calculated  for  these  models,  with 

the  parameters  (4.5),  (4.6): 

Table  4.2  Simplified  Doppler  Spread  Models 


f 

0 

k 

0 

(I).  AF=Ak 

(II).  AF=Ak  /I  +  B 

0 

20  kHz 

83.78 

6.55  Hz 

17.8  Hz 

15 

62.84 

4.91 

13.3 

10 

41.89 

3.28 

8.81 

5 

20.95 

1.64 

4.41 

, 

Clearly,  without  the  AM  contribution,  the  results  are  nowhere  near  those 
observed  empirically. 


4.2  A  Classical  Model:  Bass  and  Fuchs  et  al .  [29] 

Here  we  take  Eq.  15  of  [29]  and  use  the  same  criterion  (10.10)  for 
spectral  spread  as  above.  The  resulting  doppler  spread  for  this  "classi¬ 
cal"  case  is  accordingly  (with  Ig+p  =  l1 


a1-  +  c/  =  /  2o  here) 
x  y  x  ’ 


log  C  =  1.468, 


(4.9) 


where  as  before,  cf.  (4.1),  it  is  the  interface  motion  which  causes  the 
doppler,  and  doppler  spread.  [If,  instead,  we  consider  the  particle 
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motion  adjacent  to  the  interface  to  cause  the  doppler,  as  Bass,  Fuchs 
et  al .  [29]  do,  then  we  must  replace  a  in  (4.9)  by  unity,  e.g.. 


AFi+F  *  2ko 


/it 


log  C  . 


(4.9a) 


However,  it  is  the  moving  wave  interface,  with  its  very  large  pcQ 
contrast  vis-a-vis  that  of  the  fluid  flow  of  the  particles  producing  the 
wave  motion  ultimately  (with  pc  ~1,  where  the  net  fluid  flow  vanishes), 
which  clearly  dominates.  Moreover,  without  the  factor  aQz  here  in  (4.9), 
the  resultant  AF  is  much  too  large,  cf.  Table  4.3  f f . ] 

For  the  present  experiment,  using  again  the  PM  spectrum  (1.1)  and 
the  data  of  (4.5)  above,  we  get  directly 


AF D , r  =  0.235  k  (and  AF'  =  1.51  k  ),  (4.10) 

B+r  0  B+F  0 

shown  in  Table  4.3  and  Fig.  4.1. 

Table  4.3  AFg+F;  AFg+[- 


f 

0 

k 

0 

AF B+p ,  Eq.  (4.9) 

AF g+p ,  Eq •  (4.9a) 

20  kHz 

83.78 

19.6  Hz 

1.26  102  Hz 

15 

62.84 

14.7 

94.6 

10 

41.89 

9.80 

63.1 

5 

20.95 

4.90 

31.5 

5 .  Resul ts_  and  Conclusions 

In  the  preceding  Sections  we  have  shown  that  our  theoretical  sur¬ 
face  scattering  model,  employing  the  wind-generated  soliton  ensemble 
postulated  by  the  author  in  recent  work  [l]-[3],  accounts  very  well  for 
the  observed  data  [2].  This  is  true  not  only  for  the  measured  back- 
scatter  cross  sections  .  ,  cf.  Fig.  2.2,  but  for  the  mean  doppler 

shift,  <5f  ,  Fig.  3.1,  and  the  doppler  spread,  AF,  in  Fig.  4.1.  Indeed, 
a  noteworthy  feature  of  the  model  and  the  analysis  based  on  it  is  the 
excellent  internal  consistency  demonstrated  in  the  results  calculated 
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from  it.  For  example,  the  parameters  indicated  by  the  fit  of  theory  and 
experiment  in  the  backscatter  evaluation,  e.g.,  for  instance,  are 
consistent  with  the  soliton  ensemble  (cf.  Sec.  9)  and  the  approximations 
involved  in  determining  the  doppler  spread  component  arising  from  the 
wave  motion  of  the  large-scale  surface  [Sec.  8.2].  Similarly,  our  model 
can  account  for  the  observed  doppler  spread  AF,  which  is  noticeably  dif¬ 
ferent  from  the  theoretical  predictions  based  on  classical  approaches 
(Sec.  4.2),  and  indeed,  from  analyses  here  which  omit  the  critical  ampli¬ 
tude  modulation,  caused  by  the  random  tilting  of  the  large-scale  surface 
component. 

We  emphasize  that  the  soliton  surface  model  used  here  is  not  con¬ 
clusively  established  by  the  experimental  data  [2].  Nevertheless,  it 
is  both  plausible  physically  (cf.  Secs.  1,  9)  and  capable  of  predicting 
simultaneously  such  key  quantities  as  (back)-scatter  cross  sections 
(d^),  doppler  shift  (6f  ),  and  doppler  spread  (AF)  over  a  wide  range 
of  frequencies  (5-20  kHz  here,  and  beyond).  Indeed,  with  such  a  model 
we  have  also  demonstrated  that  the  underwater  acoustic  measurement 
techniques  of  Roderick  et  al .  [2]  are  fully  capable  of  quite  accurately 
determining  these  key  quantities  (6^),  AF).  (We  expect  similar 

results  for  the  intensity  spectrum,  as  well.)  These  results,  of 

course,  are  pertinent  to  the  particular  physical  conditions  noted 
earlier  (Sec.  1),  namely:  (i),  an  effectively  "bubble-free"  environment 
in  the  near-surface  region  [cf.  discussion  in  Sec.  2];  ( i i ) ,  high  fre¬ 
quencies  (i.e.,  large  Rayleigh  numbers);  and  (iii),  small  grazing  angles. 
When  wind-wave  interactions  are  such  as  to  produce  significant  bubble 
densities  and  distributions,  the  latter  will  dominate  and  the  soliton 
surface  mechanism  will  not  be  observable. 

It  is  instructive  to  compare  our  present  model,  vis-a-vis  experi¬ 
ment,  with  the  various  other  models  proposed  earlier  and  currently  (see 
[1],  Sec.  3.2,  for  example,  and  [5],  [9]).  We  do  this  qualitatively  in 
Table  5.1. 

The  classical  models  [(1),  (2)  in  Table  5.1]  are  almost  completely 
inadequate  to  account  for  the  observed  values  of  6fQ ,  AF  under  the 

conditions  of  the  Roderick  experiments  [2],  which  operated  in  a  "bubble- 


Table  5- 1  Comparison  of  Scatter  Model s  with 
_ (High  Frequencies  and  Small  Angles) 


Model 

^inc 

6f 

0 

AF 

Remarks:  Physical  Mechanisms 

1.  Classical:  grav.  +  capil.  waves — 
no  bubbles;  [9],  [18],  see  [l], 

Sec.  3. 

No; 

p(  15-20  db) 

too  low 

No; 

(see  Fig .  3.1) 

No; 

too  low; 

(cf.  Fig.  4.1) 

Linear  wave  generation: 

dispersive;  staurated 
spectrum,  cf.  (1.1) 

2.  Same  as  1.--C+ C(UT,  var  ,  etc.) 

( [13]- [15] ;  Sec.  1.2  here. 

Eq-  (1.1).) 

Yes 

No ; 

(Fig.  3.1) 

No; 

(Fig.  4.1) 

too  low 

C  made  larger,  cf.  Eq.  (1.1); 
dispersive;  Linear  wave 

generation;  unsat.  spectrum 

3.  Classical  +  assumed  bubble  layer; 
[5],  [6],  [9];  Sec.  2  here. 

Yes 

No. 

(no  doppler 
shift) 

No. 

AM  masked 

Not  appropriate  when  bubbles 

are  ignorable,  cf.  Sec.  2; 

dispersive;  depends  on  wind-wave 
interactions;  lin.  wave  generation 

4.  Soliton  Surface  Ensemble 

(+  classical):  [l];  this  Report. 

Yes 

(Fig.  2.2) 

Yes 

(Fig.  3.1) 

Yes 

(Fig.  4.1) 

Nonlinear  mech.  Consistent 

with  wind-wave  generation, 

[3],  [4],  [32].  Nondispersive. 

*High  frequencies,  small  grazing  angles,  etc.,  cf.  Secs.  1,  2.  ,  and  [1],  and  "bubble-free"  regimes,  cf.  last  part 
of  Sec.  2. 


ro 
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free"  regime  [cf.  remarks  at  end  of  Sec.  2  above].  In  the  unsaturated 

cases  [C-*-C(U  ,  var  U  ,  ...)  as  go  +  °°,  cf.  Sec.  1.2]  it  may  be  possible 
a  a  s  (0) 

to  achieve  the  observed  values  of  ck  ^  with  this  model,  but  the  model  is 
incapable  of  providing  the  observed  values  of  doppler  (6fQ,  AF).  A 
similar  situation  arises  with  rrtodel  (3),  in  addition  to  which  there  are 
no  observed  bubbles  effective  in  masking  the  scattering  surface,  for  the 
environmental  condition  of  the  experiments  involved.  Models  (l)-(3)  are 
all  dispersive  in  their  surface  wave  production,  which  recent  work  ([3], 

[4],  [25],  [32])  suggests  is  not  the  dominant  effect  at  these  frequencies. 

For  "bubble-free"  regimes,  our  present  model  (4),  which  involves  non- 
dispersive  surface  propagation,  on  the  other  hand,  is  entirely  able  to 
account  for  the  observed  data  and  predict  other  values,  at  different 
frequencies,  grazing  angles,  etc.  (consistent,  of  course,  with  the  high- 
frequency,  small -angle  conditions  postulated  throughout). 

To  sum  up,  we  have  followed  the  philosophy  of  constructing  a  plaus¬ 
ible  scattering  mechanism  which  is  strongly  suggested  by  the  governing 
nonlinear  physics  of  wind-wave  interactions  ([3],  [4],  [32]  and  Sec.  9  here). 
We  then  use  the  data  from  careful  experiments  [2],  with  particular  atten¬ 
tion  to  the  "ground-truth"  involved,  viz.,  geometry,  wind  speeds,  cur¬ 
rents,  etc.,  to  determine  the  numerical  value  of  the  model  parameters. 

With  these  available,  the  model  becomes  predictive  for  other  operating 
conditions,  under  the  same  general  constraints  of  small  grazing  angles, 
high  frequencies,  and  "bubble-free"  nonsurface  environments. 

Our  scatter  mechanism  is  naturally  described  by  a  three-parameter 
model  (o^,  o<.,  cs,  cf.  (9.7)  or  (9.9)),  in  the  covariance  or  spectrum, 
unlike  the  classical  cases,. which  involve  two-parameter  models  (C,  U  , 

d 

cf.  (1.1))  in  the  point-spectrum,  for  example).  This  added  versatility 
helps  the  model  to  account  for  the  variety  of  different  observed  data 
and  enhances  its  capabilities  for  prediction,  consistent  with  the  inher¬ 
ent  physical  mechanisms  which  govern  the  scattering  process.  Although 
the  soliton  surface  ensemble  postulated  in  our  current  treatments  [1],  [3], 

[4],  [32]  has  not  yet  been  directly  observed,  our  new  results  (Secs.  2-4) 
taken  in  conjunction  with  the  wind-wave  interaction  model  [4],  [8],  [32],  pro¬ 
vide  inferentially  very  strong  support  for  the  existence  of  this  surface 
scattering  mechanism.  As  can  be  seen  from  Table  5.1,  the  various  "per¬ 
turbations"  of  the  classical  linear  surface  model  fail  to  amount  ade¬ 
quately  for  the  observed  phenomena. 
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Part  II.  Analysis 


To  determine  the  role  of  doppler  produced  by  the  moving  wave  surface, 
we  must  first  obtain  the  intensity  spectrum  of  the  scattered  returns. 
Accordingly,  our  analysis  is  based  on  the  results  of  Section  7  of  [1], 
in  particular  Sec.  7.3B.  Here  we  need  to  make  a  variety  of  modifications 
to  include  the  effects  of  the  various  doppler  components  of  our  physical 
model  (cf.  Sec.  3),  specifically  to  account  now  for  the  temporal  as  well 
as  spatial  variability  ofthe  complex  moving  wind-water  interface.  As 
before  [1],  we  are  concerned  here  only  with  the  cases  of  high  frequencies 
(fQ>5  kHz)  and  small  grazing  angles  (0oj>65°),  so  that  the  "facet"  or 
geometrical  optics  terms,  along  with  the  Bragg  scatter  from  the  capillary 
components,  may  be  ignored  vis-a-vis  the  surface  soli  ton  scatter  terms 
(Secs.  1,  2  of  [1]  and  Secs.  1,  2  here). 

Although  the  correlation  distance,  ,  of  the  small-scale  scatterers 
is  small  compared  to  that  of  the  large-scale  (gravity)  wave  surface,  we 
cannot  now  neglect  the  temporal  correlation  effects  in  the  latter  in 
that  they  contribute  to  doppTer  spread,  as  we  shall  see  below.  Moreover, 
in  addition  to  the  various  random  frequency  modulations  produced  by  the 
movement  of  the  small-scale  surface  components,  we  also  cannot  neglect 
the  amplitude  modulation  of  the  Bragg  scatter  wave  surface  components 
by  the  large-scale  gravity  waves,  which  is  manifested  in  the  time- 
variability  of  the  "tilt-factor"  Although  this  latter  is 

comparatively  slow,  it  exhibits  an  important,  frequency-independent  con¬ 
tribution  to  the  doppler  spread ,  which  is  consistent  with  Roderick's 
data  [2],  Fig.  14,  and  Fig.  4.1  here.  Finally,  we  must  include  the  mean 
doppler  components,  which  represent  the  doppler  shifts  (6f  )  arising 
from  the  steady  movement  of  subsurface  and  surface  current,  and  soliton 
movement  on  the  moving  wave  surface  itself. 

With  the  above  in  mind,  we  may  apply  (7.28)  in  (7.34),  (7.34a), 

(7.35)  of  [l],  bypassing  (7.31)  and  noting  that  now  in  (7.34a) 

depends  on  the  time-difference  t  =  t^  -  t^  as  well,  as  on  vertical  angle 
0Qy.  With  the  help  of  [23]  for  the  deterministic  doppler  terms,  and 
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remembering  that  the  intensity  spectrum  is  given  by 


tfy(f)  =  2Ft{MY(x)}  =  2  J  e‘1aJTMY(i)dT  =  VY(f), 


0)  =  2TTf, 


(6.1) 


cf.  [24],  Eq.  (3.42)  et  seq.,  we  may  write  the  extended  version  of 
(7.34),  (7.34a)  here  for  the  received  incoherent  scatter  return,  X-<X>, 
as 


where 


t  =  t2  -  tv  f>0~. 


j 

-2acooc0T0  ("geometric"  or  spreading  factor,  cf. 
w  -  - 2~2 — 2~  (7.11a)  [1],  and  Fig.  5.2,  [1]; 


<4”>‘ RoTRoR 


To  *  <RoT+RoR)/co- 


Ko<T>in'/  “in‘f'>e'1“'Tdf'i  “in<f'>Ef 


Sin 


f1  =  f  -  f  ,  for  narrow  band  signals, 
cf.  (7.11b),  [1], 


and  also 


(6.2) 


(6.2a) 


(6.2b) 


u)  =  2irf  ;  f  =  "carrier"  or  central  signal  frequency, 
0  0  0 


k„  - 


2"nf  /c  =  w  /c  , 
0  0  0  0 


\  (6.2c) 


where,  as  before,  cQ  is  the  (average)  speed  of  sound  in  the  water  medium. 
(In  this  definition,  (6.1)  above,  of  spectrum,  only  nonnegative  fre¬ 
quencies  are  defined.) 

Specifically,  M^(x),  (7. 34a, b),  (7.35),  [1  ]  is  now  modified  to 

be,  for  these  cases  of  high  frequency  and  small  grazing  angles, 
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Here  we  note  the  following  changes  introduced  into  (7.35),  [l]: 


(D  exp(ik02a0*v;DT): 


© 


-A-j-Ar  /4 


This  represents  the  doppler  shift  (~  v^), 
stemming  from  the  yelocity  components  of 
various  moving  scattering  elements,  etc. 
The  anatomy  of  v^  is  described  in  Section 
7  ff . :  it  consists  of  a  mean  term  (=  dop¬ 
pler  shift)  and  a  fluctuation  component 
Ayop,  which  produces  a  random  phase  modu¬ 
lation  and  spectral  spread.  The  average 
<  >P  is  over  these  random  terms. 


(6.4a) 


This  term  is  I_(Ar),  (6.40),  [l],  and 
stems  from  the  fact  that  1^(0)  if1  (7.35), 
[l]  is  replaced  by  I^Ar),  to  prevent 
possible  singularities  in  the  integration 
over  At  in  (6.3).  [As  we  shall  see  later, 
cf.  Sec.  10,  A-j.  proves  to  be  sufficiently 
small  and  the  integrations  devoid  of  sin- 

/-V 

gularities,  that  A^  =  0  without  notice¬ 
able  effect  on  our  results.] 


(6.4b) 


®  KN(o)®r-T): 


This  replaces  i nc <— o >  ( 7 - 35 ) .  Cl], 

is  given  in  (7.28),  [1],  and  represents 
the  above-mentioned  amplitude  modulations 
by  the  large-scale  (gravity)  wave  surface. 
This  covariance  function  is  evaluated  in 
Sec.  8.3  ff. 


(6.4c) 
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„  /  2ikn^o’(S?"Ii)rA 

©  (e  /G  =  F2(-2ikoao2,  2ikQao2;  Ar,x)Q: 

=  exp  |-4a^kVz  [1  -  PG(Ar,x)]j  , 

for  the  gaussian  gravity  wave  elevation 
Cg;  (we  neglect  the  small  doppler  cor¬ 
rections  noted  in  (7.29),  [1]).  This 
is  appropriately  approximated  in  Sec. 
8.2  ff. 


(6.4d) 


(6.4e) 


Here  Kc(Ar,x)  is  the  space-time  covariance  of  the  soliton  ensemble  which 
provides  the  effective  scattering  mechanism  in  our  present  model,  in¬ 
cluding  the  (mean)  doppler  shift  due  to  Bragg  scattering,  cf.  Sec.  7  ff. 
The  explicit  structure  of  is  developed  in  Sec.  8.4  ff.,  based  on  [3], 
[4].  The  vector  quantity  2a.  is  given  by  (5.24),  [1],  and  for  the  back- 
scatter  geometry  considered  here  (cf.  Fig.  2.1  and  Figs.  2.2a,  5.1  of 
[1]),  becomes  explicitly 


2a  =  2i  cosd)  xsine  T  +  2?  sfnd)  ,  sine„x  +  2?  cose  T * 
~o  ~x  YoT  oT  ~y  voT  oT  oT 


(6.5) 


Other  pertinent  quantities  in  (6.3)  are: 


A1 

—  =  A 
2  REF 


=  reference  area,  Eq.  (2.11),  [l],  used  in  the  defi¬ 
nition  of  the  scattering  cross  section,  cf.  Sec. 
2.2  of  [1]. 


gT,  gR  =  beam  "gains,"  cf.  (6.17),  [ l] - 


(6.6a) 

(6.6b) 


The  various  approximations  on  which  (6.2),  (6.3)  above  are  based 
are  cited  in  Sec.  5.5,  [l}.  The  principal  ones  are: 

(i).  (local)  stationarity  and  homogeneity  of  (at  least)  the  large- 'A 
scale  wave  surface; 

(ii).  narrow-band  signals; 

(iii).  far-field  (i.e.,  Frauenhofer  region)  operation; 

(iv).  small  dopplers  and  slowing  moving  surfaces  (vis-a-vis  cQ); 

(v).  gaussian  statistics  for  s(r,t)G,  as  noted  in  (6.4c); 

(vi).  narrow  "beams"  (or  at  least  one  (T  or  R)  narrow  beam); 
cf.  Sec.  6,  Cl]. 


(6.7) 
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(For  others,  see  the  derivations  in  Sec.  7.3,  [1],  leading  to  (7.35).) 

Our  immediate  next  steps  are,  accordingly,  to  construct  a  doppler 
model  (Sec.  7);  to  evaluate  F^q,  (6.4d,e),  appropriate  to  the  integrand 
of  (6.3),  (Sec.  8.2);  to  obtain  (Sec.  8.3)  an  appropriate  structure 
for  ^(0) >  cf.  (6.4c);  and  then  to  construct  the  surface  soliton  covari¬ 
ance  K^,  (Sec.  9),  in  the  light  of  our  most  recent  work,  [3],  [4],  which 
supercedes  the  initial  (and  limited)  model  originally  presented  in  Sec. 

3E[  of  [1].  Finally,  in  Section  10  ff ■  we  evaluate  the  doppler  spread, 

AF,  for  our  model,  summarizing  in  Section  4.2  earlier  corresponding  results 
based  on  the  purely  "classical"  model  where  the  scattering  mechanism 
is  provided  by  the  gravity-capillary  surface  alone. 

7.  Thg  Mean  Doppler  Shift:  Analysis 

The  deterministic  velocity  v.oD  in  Eq.  (6.4a)  represents  the  sum  of 
the  mean  (or  transport)  velocities  associated  with  the  various  moving 
fluid  elements  of  the  physical  environment  of  Roderick's  experiment  [2]. 
Specifically,  these  are:  (1)  a  noticeable  deep  layer  current,  extending 
from  the  surface  to  the  bottom,  o(21  cm/sec),  and  (2)  the  surface  drift 
layer,  produced  by  the  wind  upon  the  wave  surface.  In  addition,  there 
is  a  doppler  shift  generated  by  the  sol i ton-ensembl e ,  in  our  present 
model,  which  constitutes  the  moving  water-atmosphere  interface.  This  is 
the  effective  Bragg-scatter  component  here,  which  as  we  have  shown 
earlier  (Cl],  Secs.  2,  3),  and  Fig.  2.2  of  Section  2  above)  is  the  domi¬ 
nant  contributor  vis-a-vis  the  ultra-gravity  and  capillary  elements  of 
the  wave  surface. 

The  doppler  velocities  for  the  elements  of  the  moving  mechanisms 
(1)  and  (2)  are  represented  by 


(7.1) 


which  we  can  rewrite  as 


(7.2) 


A\/q  -  j.c(vc"vc)  +  ld^Vd”Vd^  =  velocity  fluctuations. 


-kP  vd  d 
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The  deterministic  portion,  ^v  p,  is  responsible  for  doppler  shift  contri¬ 
butions,  while  the  velocity  fluctuations,  AVp,  contribute  to  the  doppler 
spread,  AF,  cf.  Sec.  8  and  Sec.  10  ff.  Figure  7.1  shows  the  present 
doppler  geometry,  along  with  that  of  the  surface  wave  number  =[J</ 2tt  = 

(v,<j>v)  =  (vx,v  )],  and  the  transmitter-receiver  configuration  for  back- 
scatter  here.  It  is  assumed  in  our  present  model  that: 

(i).  all  doppler  velocities  lie  essentially  in  the  plane  of  r,  =  0 
(equil.  wave  surface)  or  planes  parallel  to  it; 
all  dopplers  are  small:  lvn|  «c  ,  etc.; 

-tt/2 tt/2 :  waves  can  go  forward  in  the  direction  of 
the  wind,  over  an  angular  domain  it  radians  in  extent; 
the  drift  layer  velocity  is  essentially  in  the  mean  wind 


(11). 

(iii)- 


(iv). 

(v). 


direction,  i.e.,  <j>^  =  $w,  cf.  Fig.  7.1; 

the  surface  wave  numbers  jj_lie,  essentially,  in  the  (x,y)-  or 
(?=0)-plane,  cf.  Fig.  7.1. 


(7.3) 


Figure  7.1:  Geometry  of  the  fluid  transport  dopplers  (deep  current 

and  surface  drift),  backscatter,  and  surface  wave  numbers. 
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The  unit  vectors  of  (7.1),  (7.2),  cf.  Fig.  7.1,  are  given  by  the 
following  relations: 


i  =  i  cos<f>  +  i  sirKf)  ;  i,  =  i  cos$  +  i  sin<t>  , 
~c  "X  c  ~y  c  "'-d  '~x  w  w 


(7.4a) 


while  for  subensembles  of  solitons  moving  in  the  direction  ic,  we  have 


ic  =  i  cosa  +  i  sina  . 

w*x 


(7.4b) 


The  solitons  comprising  the  wave  surface  interface  all  move  in  the 

direction  of  the  wind,  i.e.,  -Tr/2<a-$  <tt/2  here,  just  as  the  gravity- 

w 

capillary  surface  waves  do,  and  like  the  latter,  they  move  ani sotropical ly 
on  the  surface,  subject  to  some  distribution  law 


w, (a)  =  A  cosa  , 
I  n 


-7T/2<a  -4  <  it/2,  A  =  <cosnct>  ^ 
'  w  n  a 


(7.4c) 


typical ly. 


7.1  Resonant  (or  Bragg)  Scatter:  General  Doppler  Shift,  I 

Equation  (6.2),  with  (6.3),  is  precisely  the  dominant  Bragg,  or 
resonant  scatter  term,  at  the  high  frequencies  and  small  grazing  angles 
of  interest  here. 

To  demonstrate  the  doppler  shift  6f  explicitly,  let  us  rewrite  (6.3) 
more  compactly  as 


/n\  rr  /  2ik  a  *v  -t\  2 i k  a  *Ar 

^RT  (T )  =  B  ^  KE(Ar,T)(e  0-0~0D)De  °~°  d(Ar)  (7.5) 


[  r] 


where 


=  K  ^o)  (^’T)F2(Ar,T)GKs(Ar,T)e 


-A^-Ar  /4 


(7.5a) 


is  the  product  covariance  of  the  various  covariances  of  the  effective 
scattering  surface  £,  cf.  (6.3),  (6.4).  This  effective  scattering  sur¬ 
face  E  is  a  "product"  surface,  not  simply  formed  as  a  linear  supposi¬ 
tion  of  surface  wave  components,  generally.  Nevertheless,  the  integral 
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in  (7.5)  is  just  the  wave  number-time  intensity  spectrum  associated  with 
Ki  -  Kz-<exp(2ikta0-voDT»[),  viz.: 


:<VT>  5 


FArl  V 


■  II 

[Ar] 


Ky(Ar,x)e 

Z_i  mvv* 


2ik  a  ■ 
o~oi 


-oDT 


(7.6) 


2ik  a 
o-oi 


2ik  a  - 
oo 


Ar 


d(Ar ) 


OO/B, 


where  we  have  used  (7.2)  to  extract  the  mean  doppler  term  (~  vQp)  from 
the  doppler  average  <  explicitly.  Here,  in  particular,  we  write 


~~ol  ^o^x,y: 


(7.7) 


from  (6.5):  there  is  no  z-component  of  a  in  v,  since  Ar  lies  solely 

w\Q 

in  the  (x,y)-plane.  In  general ,  and  AVq  may  depend  on  position. 

Next,  we  use  the  fundamental  relation  between  the  point-(intensity) 
frequency  spectrum  wl  of  an  anisotropic  surface,  and  its  directional 
covariance  function,  K£(At,t),  [26]: 

=  27ifs,  (7.8) 


LU  t  \ 

K£(Ar,x)  =  l  ^cos[iost  -  (fs)-Art)  dfg. 


wi  th 

J& ^ •  Ar  =  K  i Ar  cos (4)  -a),  (7.8a) 

(a ) 

in  polar  coordinates.*  Here  1C,  (f  )  embodies  the  dispersion  law  (re- 

_i  “  s 

lating  wavelength,  (~  K^,,)  to  wave  frequency  f  )  of  the  "product" 
surface  z1,  whose  covariance  of  elevation  (5  , )  is  given  by  (7.5a)  above 
including  the  doppler  term,  cf.  7.6.  [in  general,  ^  ,  is  complex  in  this 
formulation.  However,  without  much  loss  of  generality,  if- we  can  regard 


*Although  z  is  generated  nonlinearly,  cf.  (7.5a),  it  can  be  spectrally 
decomposed  linearly,  i.e.,  by  the  usual  Fourier  methods. 
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as  independent  of  Ar,  and  <exp(2ikQa0^*AVpT )>^  as  real;  then  , 
is  also  real-positive.] 

Accordingly,  we  shall  assume  here  and  subsequently  that  vQp  and 
Avn  are  independent  of  Ar.  Then,  applying  (7.8)  in  (7.6)  gives 


:(VT>  = 


2 i  k  i  •  A v p.T\  00 


o—ol  *~D 


i  [to  +2k  a  .  •  v  n  ji 
s  o~ol  -oDJ 


•<6(v  -F  (a,f  )/2tt)>  df  , 
~o  s  '  a  s 


(7.9) 


where  v  is  again  given  by  (7.7).  Now  Wy  is  real -positive.  [In  fact, 
it  may  be  obtained  by  the  Wiener-  Khintchine  theorem  ((3.42),  [24]) 
applied  to  (7.5a)  on  setting  Ar  =  0  therein,  viz.: 


oo 

Wz  =  2Ft(Ke>  =  2  /  Ks(0,T)K^(0)(0,x)F2(0,T)Ge"1WT  dx , 


(7.9a) 


(or  to  (7.8)  with  Ar  =  0),  where  F2_q  is  given  by  (6.4d,e).]  Now 
JFj,  =  JK^  =  (f$ ) ,  and  the  precise  form  of  jcf  ^  depends  on  the  dis¬ 

persion  law  appropriate  to  the  ("product")  wave  surface,  E,  in  question. 
Equation  (7.9)  is  consistent  with  waves  in  the  "forward"  or  wind-direction 

only,  e.g.,  -xt/2 < cx < xr/ 2 ,  cf.  Fig.  7.1  and  (7.9c).  Also,  we  have 

w 

f/_ ( f  )  =  f/y(-f  )  formally,  although  physically  only  f  >  0-  represents  real 

2j  S  2-1  S  S 

wave  frequencies. 

Next,  it  can  be  shown  [27]  that 


6(^"^s(a’fS  )/2TT)  =  6(V  "If  C°Sa)  (Vy 


=  ^  5(v-  K^ff  s  )/2tt)  6(  4>v  -  a) 


2tt i p •  (v-K  / 2tt  ) 
e  ^ 

(7.10) 


dp 


in  both  rectangular  and  polar  coordinates  [v=|vj  ,  withv  =  (v,<|>. .)]. 
In  addition,  we  have  directly  the  general  relation 


6(v  -  K ^,( f s ) / 2 xr)  -  6(f s  -  G^v)) 


Gz(v)  = 


K 


(7-11) 
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so  that  the  polar  form  of  (7.10)  becomes  directly 


6(v  -  F  /2tt)  =  - 


df 


dv 


6(f$  ■  ^(V)  * 


df 

dG 

and 

s 

- 

_ s_ 

dv 

dv 

(7.12) 


The  dispersion  law  here  is  formally  expressed  by  (7.11):  G(v)  =  f  , 
in  this  inverse  form. 

Applying  (7.12)  to  (7.9)  gives  at  once 


(7.13) 


In  addition,  we  have 


<6U 


■a)>a  = 


(cc) 


w1(a)6(<J>v 


-a)  da  = 


w.  (cfi 

1  V, 


)  ;  •  -tt/2  <  a- 


a 


$  ,<  tt/2  , 

W 


(7.14) 


where  w^(  )a  is  the  pdf  of  "wave"  directions,  as  noted  above  and  else¬ 
where,  cf.  Eq.  (7.4)  et  seq.  Our  general  result  (7.13)  shows  at  once 
the  contributions  to  the  (angular)  doppler  shift,  e.g., 


6(jo  =  2ttG  (v  )  +  2k  a.,  *v  n 
o  2  o  o^ol  ~oD 


(7.15) 


The  first  term  of  (7.15)  is  the  Bragg  scatter  component  and  the  second 
term  embodies  that  due  to  the  fluid  flow.  The  averaged  term  ,  <  >p,  in 
(7.13)  contributes  to  the  doppler  spread  [cf.  remarks  in  Sec.  7.3  ff .] . 

7.2  Doppler  Shift  II:  Various  Dispersion  Laws 

To  obtain  a  specific  doppler  shift  we  must,  of  course,  know  the 
dispersion  law  of  the  wave  surface  2.  As  an  example,  we  consider  the 
following  wave  types: 
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I.  Gravity-Capillary  Waves: 


(7.16a) 


“s  =  1 |gKs  +  f  Ks  '  or  G<v)  =  f 2?  +  2”3  f  =  V 

with  2“  the  surface  tension  parameter  =  74  dynes  cm  }  p.  176  of  [26]; 

p  =  density  of  water  (g/cirZ). 
w 

la.  Gravity  Waves  (primarily):  K 

w2  =  gK_,  K  =  w2/g;  v  =  ^  =  w2/2^g,  or  (7.16b) 

S  S  s  S  £7T  s 


fs  =  G(V)  =  v/g^/2n 


(7.16c) 


lb.  Capillary  Waves  (primarily): 


ws  =  (r/Pw)K^/2,  or  G(v)  =  W2ttv3^7pw  =  f$ 


(7. 16d) 


II.  -Soliton  Ensemble  (Non-di spersi ve) : 


f  =  c  v  =  c  K  /2t t;  G(v)  =  c  v  =  f  .  (7.16e) 

c  s  s  s  $  s 

From  (7.16)  we  readily  obtain  the  desired  inverses  f  =  G(v),  viz.: 

la.  Gravity  Waves: 

6(v  -  af2)  =  — —  6(f  -  Zv/a)  =  6(f  -  /gv/2ir)  ;  (7.17a) 

b  2^av  s  2/2TTv/g  s 

lb.  Capillary  Waves: 

6(v  -  bf2/3)  =  -^72  6(fs  '  (v/b)3/2)’  b  =  Paj/^*(2Tr)1/3;  (7.17b) 

2b  ' 

II.  Soli  ton  Ensemble: 


(7.17c) 


which  last  we  shall  use  for  our  scatter  mechanism  here,  to  determine  the 
desired  doppler  shifts  6f  (=  6coq/ 2tt ) ,  both  for  the  soliton  surface  and 
the  "classical"  wave  surface  models  I,  Ia,b  above. 
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7.3  Soliton  Surface  "Waves" 

For  (7.16),  (7.17)  to  apply  specifically  here,  i.e.,  for  us  to  use 
these  specific  dispersion  laws,  we  must,  of  course,  have  the  corres¬ 
ponding  wave  surfaces.  This  means  that  in  (7.5a)  the  wave  surface  Z 
must  refer  either  to  models  of  the  type  I,  or  II  above:  this  means,  in 
turn,  that  K^(0)»  ^-G  must  be  essentially  independent  of  Ar,  so  that 
Ks(At,t)  is  the  directional  covariance  of  the  particular  surface,  Z,  for 
which  (7.16),  (7.17)  hold. 

For  the  present  soliton  case  we  Use  (7.17c)  in  (7.13),  along  with 
(7.6),  (7.8a),  to  write  directly 


,(0) 


"rt'MV  *  ■  D„H4>vvfo 


2-VlI/co> 


(iio  [c  |2oi  |  /  c  +2a  .  •  v~/ c  ]x} 
oL  s 1  0]_  1  o  '«-0l  ~-oD  o 


(7.18) 

with  G(vq)  =  c  v  .  [Here  D0(t[v0)  is  given  explicitly  by  all  factors  in 
(7.13)  not  containing  the  doppler  shift  term.]  The  doppler  shift  con¬ 
tributions  are  seen  at  once  to  be  the  coefficient  of  x  in  (the  exponent 
of)  (7.18),  as  we  shall  show  in  more  detail  below.  The  quantity  Dq  is 
real  and  positive,  vanishing  according  to  (7.14)  outside  the  allowed 
interval  for  a-$w:  this  corresponds  to  soliton  "waves"  moving  downwind 
only,  and  not  against  the  wind  direction. 

Our  next  step  is  to  apply  (7.18)  for  M^(x)  to  (6.2),  observing 

HI  /Q\ 

that  except  for  the  exponential  term  in  (7.18)  M^y(x!fo)  is  even  in  x. 

Since  K  (x).  =  K  (-x).  ,  cf.  (6.26),  we  can  write  the  intensity  spec- 

oinoin 

trum  (6.3)  of  the  (incoherent)  scatter  return  as 


W , 


X-<X> 


(f) 


k2G^ 


K  (x).  D  (x|<J>  ,v  lcos(w-b  )x 
o'  'in  o'  v0  o  o' 


(7.19) 


where  now 


/\ 


=  U0(l 


cs/co 


co  +  6co„ . 

0  0 


(7.20) 


[We  have  chosen  the  signs  in  the  frequency  dependent  terms  in  (7.20)  to 
conform  to  the  geometry  of  Fig.  7.1:  the  scattering  region  is  downwind 
from  the  receiver.  Thus,  for  downwi nd  here,  we  have  0<4>w<tt  ;  for 
upwind ,  we  have  tt<$w<2tt.  From  (7.20)  it  is  immediately  evident  that 
the  desired  (general)  doppler  shift  6f  is 
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6u> 


5fn  =  9 
o  2v 


o  


1 2a  .  I  c  f  /c  -  2a  i  •  v  ->f  /c  . 
~oJJ  s  o  o  ~ol  ~oD  o  o 


(7.21) 


Applying  (7.4a)  to  (7.2)  and  using  (7.7)  gives  explicitly,  for  the 
present  backscatter  geometry  (<j>  y  =  tt/2): 


(7.22a) 

(7.22b) 


which  is  the  relation  needed  in  Sec.  3,  earlier.  Note  that  for  the 

deterministic  components  of  the  moving  fluid,  i.e.,  the  deep  current 

and  drift  layer,  the  azimuthal  angles  of  observation  (d>  ,  $„)  appear 

c  w 

explicitly,  while  for  the  random  moving  i nterface ,  determined  by  the 
soliton  ensemble,  it  is  the  magnitude  of  the  solitons'  speed  (on  the 
drift  layer)  which  is  important,  without  azimuthal  angles.  [The  azi¬ 
muthal  angle  of  viewing,  <j>  here,  influences  the  observed  spectral 
intensity,  through  (4>v)  in  Dq,  cf.  (7.18),  but  not  the  doppler  shift 
component. ] 

In  addition,  from  (7.18)  and  (7.19)  it  is  at  once  evident  that 
when  AVp  =  0,  here,  so  that  Dq  is  independent  of  x,  then  Dq  is  propor¬ 
tional  to  the  "spectral  weight"  or  line  intensity  of  the  Bragg  scatter 
component  of  the  moving  surface,  selected  by  our  choice  of  illuminating 
frequency  f  .  Because  of  the  planar  character  of  the  dispersion 
relations  (7.18),  (7  17)  for  these  small -amp! i tude  surface  "waves,"  the 
resulting  frequency  spectrum  is  a  line  spectrum.  Thus,  letting  the  il- 
luminating  waveform  be  purely  sinusoidal,  so  that  K  (t).  =  (AQ/2), 

cf.  (6.2b),  we  see  that  (7.19)  becomes 


vx-<»(f) 


cw 


-  D0(WkoG(l)  T  >  cos(d3-bo)x  d- 


A 

-  ko  f  G(1)D„(»1,.v„)«(f-[V6fJ);  v„  =  |2on,  |f„/c„. 


0  TV  O' 


0  0- 


'Ol  1  0  0  ' 


(7.23) 
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as  remarked  above.  For  (narrow  band)  but  not  purely  sinusoidal  illumina¬ 
tions  (7.19)  reduces  directly  to 

Wf-[fo+5V>  =  koG<1)[lo(Vvo>  win[f-(V6fo>]’  <7-23a> 

which  has  the  spectral  shape  of  the  narrow-band  injected  signal,  cf.  Fig. 

11,  [2]. 

Actually,  as  we  shall  see  explicitly  later  [cf.  Secs.  8,  10],  the  line 
spectrum  (7.23)  is  somewhat  broadened,  because  of  the  "doppler  smear" 
produced  by  the  random  variability  of  the  various  fluid  and  interface 
velocities,  e.g.,  Av^  f  0,  (7.2),  etc.,  and  the  random  "tilting"  of  the 
large-scale  surface,  embodied  in  K^(q),  cf.  (2),  Eq.  (6.3).  The  result¬ 
ing  doppler  spread,  AF,  accompanying  the  doppler  shift  <5f  ,  cf.  (7.22) 
here,  may  be  calculated  from  Do(4>v,^0)  by  a  perturbation,  allowing  •> 
vq  +  Avq  (~  AF/cQ)  therein,  or  by  a  direct  calculation  of  (7.5)  in  (6.2). 

In  either  case  we  need  the  specific  structure  of  the  covariance  func¬ 
tions  involved,  viz.  (7.5a).  These  we  obtain  in  Section  9  ff. 


7.4.  The  "Classical"  Wave  Surface 

Using  (7.17a)  in  (7.12)  and  applying  it  to  (7.13)  gives  directly 
the  wave  number-time  intensity  spectrum  for  the  purely  gravity-wave  sur¬ 
face  (X  >2  cm).  We  have 
o 


(v  r )  =  —  .1— ^ 

Gv  o  'grav. 


4  |2ttv^  Fl(21koao1-DT,D  fc'G(‘/gv0/2n)<6(4,v-a)> 

.  i(2ir/gv0/2u  +  2^2o‘70D)t 


a 


(7.24) 


where  F^_^  is  the  average  <  >^,  cf.  (7.13a).  Using  (7.14)  we  readily 
obtain  the  counterpart  of  (7.18),  viz.: 


M^Mf  )  =  B$v  =  D  (tU  ;v  )  exp 

RT  v  1  o'grav.  E-grav.  ov  |yv  o'grav.  v 


ro)„ 


w0 

- 

r9 

0 

- , 

— 1 
O 

3{ 

C\J 

(7.25) 


where  now  the  weighting  Dq  is  explicitly 
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D  (x  I  cf)  ;v  ) 
o  |Tv  ~o'grav. 


B 

4 


II-9- 

Pvo 


Fl^2iko-oi'-iDT^D 


J7E(/gvo/2TT)Wi(<j)v)  (7.25a) 


analogous  to  (7.18). 

The  desired  doppler  shift  is  obtained  as  before,  cf.  (7.19)-(7. 22 ) , 
and  is  at  once  seen  to  be,  generally, 


Sf0  =  -[  >'9|2Vl/2,lX0  + 


(7.26a) 


and  specifically  here  for  backscatter  (with  <£  y  =  tt/2) 


Sf 


g  sine  T  2sin0 


oT  ,  oT  r- 


ttX, 


X, 


[v  sin4  +  v  .  sin$  , 
L  c  rc  d  wJ1 


(7.26b) 


including  the  surface  drift. 

For  the  complete  gravity-capillary  wave  surface  we  use  (7.16a)  in 
(7.12)  to  get 


6(v-<F/27r) 


g+3K2e 

2H+BKo 


gvQ/ 2tt  +Bvo2Tf]5(<f)v-a) ; 


3  =  e7/p  ;  K  =  2ttv  . 
'wo  0 

(7.27) 


Accordingly  we  replace  the  first  term  of  (7.26b)  by 


g  sineoT 


si  n20 


ttX 


+  1  6tt  3 


oT 


0 

for  the  complete  effect  on  6f  now.  However,  for  our  present  applica¬ 
tions,  (7.26b)  is  quite  sufficient,  since  X  >7.5  cm  here,  and  the 

o  i 

capillary  regime  does  not  dominate  until  about  A  <  (^  -  1)  cm,  cf. 
p.  178  of  [26]. 


(7.28) 


8.  Doppler  Spread,  I: 

To  obtain  the  doppler  spread, AF,  we  need  the  explicit  structure  of 
(D  >  (2),  and  (D  in  (6.3).  [It  turns  out  here  that  we  do  not  need  the 
explicit  covariance,  Kc(Ar,x),  of  the  soliton  ensemble  for  this  purpose. 

j 
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Of  course,  K<.  is  required  if  we  desire  the  (intensity)  spectrum  of  the 
scatter  return  itself,  and  it  is  needed  for  the  backscatter  intensity, 
cf.  Sec.  2.  See  Eqs.  8.9  ff.]  Accordingly,  we  begin  with: 


8.1  "Flow"  Spread 

Here  we  need  to  evaluate 


Fl^2lko-or~DT^D 


i 2k  a  ,‘AVnT 
0-0 1  " 


X 


(8.1) 


for  the  fluctuations,  Avn,  in  the  deep  current  and  surface  drift  veloci- 
ties,  cf.  (7.1),  (7.2).  Making  the  usual  (reasonable)  gaussian  assump¬ 
tion  for  each  independent  component  of  av^,  writing 


A  v 


c 


(2aoi-jc)Avc; 


Av 


(2a 


oi 


ld)AV 


(8.2) 


gives  directly 

=  ^VaV°o)  (JV'V'o 


(8.3) 


2  , 


Here  the  variance  is  specifically 


°B  £  (2“oi'lc/co)20c  +  (2SoI-V0o>V  °c  5  Avg;  °d  5  Avd  '  <8-4> 

These  random  fluctuations  produce,  as  expected,  cf.  (8.3),  a  random 

phase  modulation,  which  depends  on  the  input  signal  frequency  f  . 

In  our  present  model  of  the  soliton  surface  ensemble  we  postulate 

that  the  variance  in  the  soliton  surface  speed  c^  is  essentially  negli- 
?  2  2  ^ 

gible,  so  that  Ac^  =  0,  e.g.,  =  Ac^  *  0:  there  is  negligible  doppler 

spreading  due  to  variations  in  the  speed  of  travel  of  the  numbers  of 


the  soliton  ensemble  generated  by  "wind-wave-wind"  interactions  on  the 
large-scale  surface  (as  described  in  [3]  and  [4]).  The  major  variations 
appear  in  the  size  of  the  individual  solitons  which  comprise  the 
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hydraulic  impulses  produced  on  the  large-scale  surface,  cf.  Fig.  3.1(d), 
[4],  for  example.  [See  Sec.  9  for  additional  remarks.] 


8.2  Spreading  Due  to  the  Large-Scale  Wave  Surface  Motion, 

I.  Random  Phase  Modulation 

Next,  we  consider  the  effect  of  the  (primarily  here)  temporal  vari¬ 
ations  in  the  large-scale  surface,  primarily  the  gravity-wave  component 
(G),  as  exhibited  by  (4),  namely  F^g,  ^ •  (6.4d,e).  These  also  produce 
a  random  phase  modulation,  dependent  on  the  illuminating  frequency  (f  ). 

Using  (6.4e)  and  the  basic  relation  (7.8)  allows  us  to  write  for 
the  directional  covariance  of  the  large-scale  wave  surface  (G)  here 


2 

7gpg' 


on  pr(Ar  ,t  ) 


=  /  ^(fs)<cos( 


to  t  -  Kr •  Ar )>  df  ,  t  =  t0-tn,  Op  =  C 
S  ^  a  S  L  I  b  L 


(8.5) 


with  J<q*At  =  KgAr  cos(A<t>-a),  (Ar  =  r^-r^,  Ar  =  l^fh  A<J>  =  <f> 2"<J>i>  etc-)> 
and  Kg  (=K$)  of  (7.16),  in  the  usual  way.  Here  Vg(f  )  is  The  "point"- 
intensity  spectrum  of  the  wave  surface  elevation  c(r,t).  [The  semi- 
empirical  model  we  use  here  explicitly  for  is  the  well-known  Pierson- 
Moskowitz  result  [15],  Eq.  (1.1)  above.] 

For  the  large  Rayleigh  numbers  involved  in  the  present  study  we 
note  that  the  principal  contribution  to  the  exponent  of  (6.4e)  arises 
for  small  _Ar  and  x.  Consequently,  we  seek  a  Taylor's  expansion  of  Pq 
about ^Ar  =  0,  t  =  0,  with  as  yet  unspecified.  For  a  pdf  of  a  of  the 
form  (7.4c)  for  the  different  wave  directions,  we  readily  find  that  the 
only  non-vanishing  terms  in  the  expansion  (through  the  second  order)  are 
such  that 


ar.Pr,(Ar,t)  =  1  -  ^j-j  a2(Ax)2  +  a2(Ay)2  + 


2  2 
x  ot 


3  3, 


GKG 

where  specifically 


Or 


2  '  tAxoxx  r  +  °(Ax  ’t  ) » 


(8.6) 


=  /  k; 


:os2oe^ 


wr  df 


G  Vsin2ct7  G  s 


(>  o)-, 


(8.7a) 


2 


0 


XT 


/  u^Kg  cosa  ffg  dfs  (>  0) ; 


(8.7b) 
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4  =  4  “s  w&  dfs  ■  ?2  <> 

Applying  (8.6)  to  (6.4e)  allows  us  to  write  now 


(8.7c) 


F2-G  e 


'(2aGkoaoz)  ^ 1 " PG ^ 


-  l-(2a0zk0)2[a2Ax2+a2Ay2+T2a|]+I(2a0zk0)2xAxa 


=  e 


2 

XT 


(8.8) 


At  this  stage  we  use  the  P-M  spectrum  (1.1)  and  make  a  variety  of 

approximations,  to  simplify  (8.8)  vis-a-vis  the  small-scale  (soliton) 

2 

surface  component.  Let  us  consider  first  a  ,  (8.7b).  We  have,  since 

r\  XT 

kg  =  V9,  from  (7-16b)’ 


XT 


PM 


=  /  0)1-2- 
n  s\  g 


cosa 


r2  - 
Cg  to 


-ag 


s  a 


df 


[a  =  0.74,  C  = 


8.10-10"3] 

(8.9a) 


Oleosa 
2  TT 


4 

00  _2  -B/x4 
/  x  e  dx 
0 


Cg  cosa 
"  8*B1/4 


oo 

/  y‘3/4e-ydy; 
0 


4  — 4 

39  /V 


PM 


CU  T( 1/4 )  cosa 

a 


8Tra 


1/4 


(8.9b) 


Next,  we  use  r( 1/4)  =  4T(5/4)  =  4(0.906)  =  3.624;  a1/4  =  (0.74)1/4  = 

2  2  2 

0.93;  and  a  cos  a  law,  viz.  w^(a)  =  —  cos  a  from  (7.4c) ,  -tt/2  <a<TT/2, 

to  get  cosa  =  3/ it,  so  that  (8.9b)  becomes,  at  U  =10  m/sec, 

a 


PM 


(8.10  10“3) 

8  TT 


10(3.624)  _ 
0.93 


1.20  •  10"2  (msec)-1 


U  =  10  m/sec 
a 

(8.9c) 


Also,  from  (2.5),  the  Cox  and  Munk's  relation  [17],  we  obtain  (as  a 

2 

result  of  the  w^(a)  =  (2/tt)cos  a  assumption) 

U  =  10  m/sec. 
a 


(8.10) 
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In  addition,  for  the  strictest  case  here  at  f  =  20  kHz,  we  have 
kQ  =  83.78  rad  m  ,  2aQZ  =  2  cose  y  =  2  cos  81°  =  0.313,  so  that 

(2aozk0)2  =  [(0.313)  •  83. 78]2  =  [26. 2l]2  =  6.87  •  102  (m"2)- 


(8.11) 


At  this  point  we  use  the  (normalized)  space-  and  time-covariances 
of  the  soliton  ensemble  to  obtain  their  correlation  lengths  and  times, 
from  (9.7),  (9.13),  namely. 

2Ar  2c  x 

X-1K1^X-1^  =  e  *  X-1  ah  ’  x-l  =  ~~c^  •  (8.12) 

From  Watson,  p.  737  [30],  we  see  that  x_^  =  1.65,  so  that  using  =  5.0 
cm,  from  the  backscatter  cross-sectional  and  doppler  shift  data  of 
Roderick,  cf.  Fig.  2.2,  3.1,  appropriate  here,  we  get  for  the  correla¬ 
tion  lengths  and  times  of  the  surface-riding  soliton  ensemble 

Arc  =  PP  °h  =  4-13  ‘  10"2  m;  Tc  =  PP  Vcs  =  °'28  sec’  (8.13) 

where  we  have  used  c  =  15  cm/sec,  in  line  with  Table  3.1. 

s  2 

There  remains  the  calculation  of  a*,  (8.7c),  for  the  PM  spectrum. 

VJe  get 


,  U  =  10  m/sec.  (8.14) 

a 


Let  us  now  use  (8.13),  Ar  =  |/ax2+AY2  =  Arc,  (8.10)  and  (8.11)  to 
compute  the  purely  spatial  part  of  the  exponent  of  (8.8).  We  get 

+  <  i(2aozko)2axArc  =  ^(6-87  -102)  (5.42)10'2(4. 13  10"2)2 

=  3.17  •  10 


a?  =  C  /ir  U  2/8a1/2  =  2.09  •  10"3  U  2  =  0.209 

Q  cl  /  d 


Similarly,  we  use  t  +  t.,  Ax^-Arc,  (8.9c),  (8.11)  to  get 


(8.15) 
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i(2aozko)2xAXaxt  <  !(2aozko)2xArc°XT  =  1(6>87  102)(4.13-  10"2)(1.20-10"2)t 

,  ‘  (8.16) 


0.170  x 


and 


|-(2aozko)2a?T2  =  |-(6.87*  102)(0.209)t2 


2 

7 1 . 8t 


(8.17) 


Fort  =  o(t  )  =  0.3  secs,  cf.  (8.13),  clearly 


0. 170t 

c 


0.051 <<  71. 8(0. 3)2 


6.46 


(8.18) 


2 

so  that  we  can  very  safely  neglect  the  last  term  (~  a  )  in  the  exponent 
of  (8.8).  In  addition,  from  (8.15)  in  (8.8)  we  see  that  the  spatial 
contribution  is  negligible,  so  that  we  obtain  the  much  simpler,  approxi¬ 
mate  form  for 


-tu^T2R./2 

O 

F2-G  4  e  ; 

R:  i  (2a  a*/c  ) 

X  oz  i;  o' 

(8.19) 


The  result  is  thus  purely  a  doppler  spread  due  to  the  random  phase  modu¬ 
lation  of  the  incident  signal  (~  fQ),  produced  by  the  movement  of  the 
local  large-scale  surface  (G),  as  we  might  expect.  For  frequencies  less 
than  f  =  20  kHz,  the  above  approximations  are  even  better. 

8.3  Spreading  Due  to  the  Local  "Tilting"  of  the  Small-Scale 

Surface  Component,  II.  Random  Amplitude  Modulation 

The  third  element  producing  doppler  spread  is  the  random  amp! i tude 
modulation  generated  by  the  time-varying  slope  of  the  large-scale  sur¬ 
face  ,  on  which  the  small-scale,  soliton  ensemble  rides.  Because  of 
the  assumed  local  homogeneity  (and  stationari ty)  of  the  large-scale 
surface,  4(jr,t)g,  the  statistics  of  these  slope  variations  are  also 
essentially  independent  of  position  for  the  small  Ar  =  Arc  regions 
involved  vis-a-vis  the  soliton  ensemble. 
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The  evaluation  of  K  (g)(Ar,T)  in  (6.3)  proceeds  from  Eq.  (7.28), 

j\|  '  '  VVSA 

[ 1] ,  where  we  have  explicitly 


(8.20) 


where,  as  before  (cf.  (2.2a,b),  [  l] ) , 


4 >/(1  +  A  +  cy>’S  <=  "z  E  (I  +  4  *  ; 


nP  =  (i  c  +  i  r,  -  i 
-G  -xx  ~y  y 


.*•  n  =  (i  c  +  i  4  -  i  )n  . 

~  '—x  x  —y  y  — z  z 

2  2  -2 

For  the  small  slopes  encountered_here  (a  ~  o  -5*10  ,  cf.  (8.10)),  we 

22  *  y 

can  drop  the  denominators  (c  -»-c  ,  etc.)«l,  on  the  average,  and  so 

X  X 

write  (8.20)  approximately  as 

,  2 


(8.20a) 

(8.20b) 


K 


(0)  <(ai‘2Vnzl)  ^2’2“o/nz2^2>G* 

II 

Next,  we  use  the  following  reasonable  conditions  and  assumptions 


(8.21) 


further 

to  simplify  (8.21) ,  viz. : 

(i). 

the  pdf's  of  are  symmetrical  in  -co<x,y<0°:  all 

moments  vanish; 

odd 

(8.22a) 

(11). 

Cx  "  Cy  =  °’  sxlCy2  “  ~r,x2Cyl  =  <^CxlCyl^?x2i’y2^>  = 

0. 

(8.22b) 

(Here  1 

refers  to  t^,  2  to  t„,  =  0  on  the  soliton  scale.) 

Expandi ng 

(8.21),  subject  to  (8.22),  and  retaining  (even)  terms  of  order  no  higher 

than  the  quadratic  (because  of  the  small-slope  condition,  (8.10),  etc.), 

we  get,  for  our  present  geometry  (<f>  T  =  tt/2),  cf.  (2.4)  and  Fig.  2.1 

(and  the  assumption  of  isotropy:  a  =  a  ,  etc.) 

x  y 


,(0) 


4  /p 
oT  u 


16  {  4=)  sin\T  +  sl'n20oTcos2eoT  +  cos40otIg* 

(8.23) 
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2  2 

a  -  a 

x  y 


a  .  Under  the  further  reasonable  assumption  that  t,  ,  x, 
x  y  x  y 


are  normally  distributed  (since  is  normal),  we  have 

f t  1  =  /  C*  =  (?PUi 

x=y 


pSa)  -  44 z?f  -  (2P<ir  -  D-. 


x=y 


px=y  “  SxlCx2  • 


(8.23a) 


Then,  K  ^g^,  (8.23),  becomes 


KN(°>(~’T>  =  kn(0)(t>  4 


4  H)pJsin40oT+2(24l>}a2s1n2eoTcos20oT+cos4eoT}, 

(8.24) 

which  reduces,  as  expected,  to  (2.4)  when  x=  0,  p^,(0)  =  1. 

x-y  1 

Our  next  step  is  to  consider  the  expansion  of  p  _  (t)  about  x  =  0, 

(1)  x-y 

since  it  is  only  for  small  values  of  x  that  P^=y  makes  a  significant  con¬ 
tribution.  We  have  from  (8.5) 


(  1  \  O  p  p 

Px=y(Ar>x)g  =“/  Kg  <(cos  a~  sin  a)cos (cosx  -  Kg*Ar )> ^  i7g  dfg 


(8.25a) 


x=y 

so  that 


pl=y(0,0)G  =  -(cos2a =  sin2a)  /  w^K2  ^G(fs)dfs  (<  °°:  physically),  (8.25b) 


(>  0) 
(8.26) 


This  defines  r„,  namely, 


rG  = 


Pxly(0)G 

2! 


/  ws  ^G(fs)dfs  (>  0) 
2g  0 


(8.26a) 


when  (7.16b)  is  used. 
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If  we  use  the  P-M  model  (1.1)  for  in  (8.35b),  we  get,  with  the 
help  of  the  dispersion  relation  (7.16b), 


M) 

px=y 


,  — 9-  -  -agVu/ 

G:  PM  C  cos  a  a)$e  dfs 


(8.27) 


which  is,  of  course,  not  physically  valid,  since  at  sufficiently  high 
frequencies  (to^)  the  wave  spectrum  must  vanish,  e.g., 


5<n 

Hx=y 


G:  PM 


4/  4rr“4 

O  -ag  /a)  U 

=  -  C  cos  a  /  K  03  e  a  df  > 

os  s 


(8.27a) 


(Rather  than  attempt  to  evaluate  (8.27a),  since  w,  is  uncertain,  we  shall 
determine  I^,  or  R^,  (8.29b),  instead,  from  the  experimental  measurements 
of  doppler  spread,  cf.  Sec.  4.1. 

Using  (8.26)  in  (8.24)  allows  us  to  write* 

2 

K  (q)(t)  "  16{(3-4T2rG)o4sin40oT  +  2(3  -  ^  rG)c2sin20oTcos20oT  +  cos40oJ}, 


(8.28a) 


which  can  be  put  in  the  form 


l  -  r„  V 


fjA-iAoT  +  2oj;si  n20oTcos20oT 


G2^4.4  ??  2  4 

3axsin  0oT  +  6oxsin  0oTcos  0oT  +  cos  0oT)J 


(8.28b) 


which,  in  turn,  for  small  x  can  be  rewritten  as 


"  (0)(T) 


*qT*”/2 


*  N«)(a  te'^. 

GS  -o' 


(8.29a) 


Here,  specifically, 


RG  -  rG 


8oxs1n  9oT  +  2oxs1n  9qTC0S  9qT 

L3a4si  n40OT  +  602si  n20oTcos20oT  +  cos4e  T. 


(8.29b) 


*Since  we  are  interested  only  in  small  x,  we  write  p^2  =  (1  +  t^Tq+  )2 

.2  4  ^  7 

=  1  +  2x  rG,  dropping  terms  q(t  )• 
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with  Tg  determined  from  (8.26a),  or  (8.27a),  when  possible.  In  any  case, 
Kg  can  be  found  (as  is  done  here)  from  experiment,  by  measuring  the 
doppler  spread  at  several  different  frequencies. 

The  important  feature  to  note  here,  in  (8.29a),  is  that  the  random 
“tilt"  of  the  large-scale  surface  produces  an  amplitude  modulation  of  the 
scattered  signal,  which  is  independent  of  signal  frequency.  It  is  this 
effect  which  accounts  for  (i)  the  much  larger  observed  doppler  spread, 

AF,  and  (ii)  the  essentially  non-vanishing  character  of  AF  as  signal 
frequency  (f  )  is  decreased,  cf.  Fig.  4.1. 


9 .  A  Concise  Summary  o f_  M i d  d  1  eto  n_'  s_  Sol  i toy 
S u r f a c e -  C ornpon e n t_ Mo de  1_  [3 ] [4j 


As  noted  in  Sec.  1  above,  after  the  nonlinear  wind-wave  interactions 
have  reached  a  point  of  sufficient  development,  the  resulting  surface 
distortions,  riding  on  a  thin  wind-drift  layer,  which  in  turn  lies  upon 
the  large-scale  gravity-capillary  component,  represent  ensembles  of 
solitons,  of  varying  magnitude  and  duration.  These  surface  distortions 
are  thus  decomposible  into  a  set  of  hydraulic  bumps  (solitons),  as  observed 
by  Peterson  [3l]  and  others  [3];  [4],  Fig.  3.1,  for  example.  The  thin, 
wind-drift  layer  of  effective  thickness  h  ^  provides  the  necessary 
channel  for  the  soliton  to  be  created  and  move  upon  These  solitons  are 
then  produced  as  the  limiting  solutions  of  the  famous  Korteweg-deVries 
(KdV)  equation,  expressed  in  normalized  coordinates  by  a  relation  of  the 
form 


V  -  6??t  +?xVx'  =0'  ■  (9-D 

where  4  is  the  surface  displacement,  cf.  [10],  pp.  463-468  and  [11], 

Secs.  1.2,  5.3.  Here  the  various  solitons,  or  solitary,  one-sided 
hydraulic  "bumps,"  are  propagated  under  conditions  where  the  dispersive 
and  nonlinear  forces  in  the  (thin  layer)  of  the  medium  are  essentially 
balanced,  [10],  [11],  and  dissipative  mechanisms  (e.g.,  viscosity)  are 
very  weak  and  ignorable.  (Ultimately,  the  effects  of  surface  tension  must 
be  included  also.  These  appear  as  a  negative  component  of  the  dispersive 
forces  and  reduce  the  coefficient  of  the  dispersive  term  (~4X'X'X')  1n  the 
KdV  equation.  The  effective  drift  layer  is  now  the  soliton  surface.*) 

*D.  Middleton,  "A  Proposed  Soliton  Mechanism  in  Wind-Wave  Surface  Genera¬ 
tion  and  Scattering,"  NUSC  Tech.  Doc.,  in  prep.,  June,  1986. 
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The  (denormal i zed)  solution  of  (9.1)  can  be  shown  to  have  the  form 


C(r, t)  =  3vq  sech 


£  f>r\-(c  +vj(t-t,)] 


(9.2) 


at  some  point  where  we  are  considering  the  typical  member  of  the  ensemble, 

initiated  at  (r.,t.).  Here  v  is  the  excess  speed  due  to  nonlinear  effects, 

<1  J  „  _2 _ 

and  c  =c  +v  ,  where  c  =  /gh  rr.  Thus,  c  is  the  constant  speed  of 
s  o  o  o  ett  s  - 

the  soliton:  solitons  travel  at  constant  speed  and  their  propagation  is  non- 
dispersive,  i.e.,  c$v  =  f  ,  cf.  (7.17c),  unlike  gravity  and  capillary  waves, 
which  are,  of  course,  dispersive,  cf.  (7.17a,b).  Here  vQ  =  vqj  ,  which 
depends  on  the  decomposition  of  the  surface  into  particular  soliton  members, 
each  traveling  at  excess  speed,  vQ j ,  while  0=  h^2  /g/6.  It  can  be  shown 
that  (9.2)  can  be  wel 1 -approximated  by  [3] 


Z(r-ryt-t-)  =  ( v  Q )  j  e 


■4|r-rres<t-tJ)|2/l);j 


L_,  =  ^ 


oj 


>1/ 4.5/ 4 
■g  neff; 


oj 


h 


;oj  =  2 


eff 


v  .  (9.3) 

9  oj  ' 


Here  L  .  is  the  "wavelength"  of  the  typical  soliton,  and  may  be  regarded 

^  vJ 

as  a  random  quantity  over  the  soliton  ensemble  { j } ,  in  consequence  pri¬ 
marily  of  the  composite  random  surface  displacements  or  distortions,  riding 
on  the  drift  layer  and  large-scale  surface  component  (G):  many  different 
sizes  (~v0j)  of  soliton  are  required,  to  constitute  the  surface  distortion 
at  r.  Since  ^  is  a  limiting  form  of  solution  of  the  KdV  equation  (9.1), 
a  condition  for  (one-sided)  solitons  to  exist  is  that  (for  each  j) 


C  /h2ff  >  8,  with  L  =  A  / 2;  A  =  "wavelength"  between  e"^  points, 
o  o/  err  ooo  (g_3a) 

Otherwise  (two-sided)  cnoidal  waves  are  generated  ([10],  Eq.  100,  p.  466). 

The  resulting  soliton  ensemble  defining  the  incremental  (i.e.,  small-scale) 
surface  at  r  is  then  given  by 


4s(r,t)j  =  l  c(r-r  •  ,t-t  • )  =  \ 
j  J  J  j 


;4|r-,ycs<t-V|2/Loj 


’oj 


The  soliton  ensemble  is  then 


,  over  all  allowed  values  of  J. 


(9.4) 
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pdf 


A  reasonable  distribution  of  values  of  the  L  .  is  provided  by  the 


2  _  ,  2  .  r.  ~  _ 

ah  -  L0  ’  " 


0  =  1 


v 

TT  h 


(9.5) 


Because  of  the  many  forces  influencing  the  local  surface  at  r  a  Central 
Limit  Theorem  argument  may  be  invoked  to  suggest  the  hal f-gaussian  pdf 
(9.5)  (since  Lo>0).  Moreover,  this  pdf  favors  solitons  with  lengths 
shorter  than  (L^)^  (=  a^) ,  although  comparatively  large  values  of  LQ 
(>  )  are  not  improbable,  and,  indeed,  are  needed,  to  provide  sig¬ 

nificant  backscatter  returns  at  the  longer  Bragg-scatter  wavelengths 
(Xq/2  ~15  cm,  at  f  =  5  kHz;  Xq/2 = 7  5  cm,  at  fQ  =  10  kHz),  cf.  Fig.  2.2. 
[Another,  supporting  feature  of  our  choice  of  pdf  (9.5)  is  the  perhaps 
fortuitous  fact  that  it  leads  to  a  wave  number  intensity  spectrum,  W9(k,0)c, 

-A  C.  ~  O 

which  is  o(k  )  as  k  becomes  large,  cf.  Eq.  9.15  ff.J 

An  important  characteristic  of  solitons  is  the  unusual  feature  that 
they  remain  undistorted  when  interacting  with  one  another  [l 1 ] .  This 
suggests  that  the  soliton  ensemble,  a  representation  (J)  of  which  is 
given  by  (9.4),  may  be  described  by  a  poisson  process  in  time  (at  any 
point  r):  superposition  does  not  change  the  individual  "wave"  forms. 
Accordingly,  we  can  show  quite  readily  [3]  that  the  temporal  covariance 
associated  with  9.3  is  given  by 

2  2 
“24y /Lq 

KS(0,t)  =  °s  <L0>  ’  Ay  =  CST’  ^9‘6^ 

o 

where 

ks(o,o)  =  4  =  (9-6a) 


is  the  resulting  mean  intensity  of  (all)  the  soliton  components  at  (r4,t.) 

J  \J 

[(r,.,t.),  does  not  appear  in  (9.6)  because  of  our  assumption  of  stationarity 

vJ  J  kJ 

and  homogeneity];  7T7  is  the  average  number  of  solitons  per  unit  length, 
traveling  in  the  direction  of  the  wind,  i.e.,  downwind  generally.  Applying 
(9.5)  gives,  for  the  point-covariance  (9.6), 
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Ks(0, t)  =  a  5 


2c$  T 


2V 


2  _  |  2 
o,  -  L 
h  o 


(9.7) 


(where  is  a  modified  Bessel  function  of  the  second  kind;  note  that 
xio  xK}(x)  =  !•)•  Here  cs~co,  since  c0>>v0>  cf.  remarks  (1.1),  p.  54  ff. 

The  associated  point-intensity  spectrum  f7^(fs)  is  now  readily  ob¬ 
tained  from  the  Weiner-Khintchine  theorem  ([24],  p.  143): 

00  -id)  T 

J7s(fs)  =  2  /  Ks(0,T)e  dT  (9.8a) 

-oo 

2  °° 

2  a  1  o  ^  r 

=  — -  J  xK  (x)cos  ax  dx;  a  =  d  o, /2c  (9.8b) 

so1  .  s  n 


(t\o,  o2 

h  S 


1- 


Ksoh\2l3/2 


(9.9) 


since  w/c  =  Kc,  (7.17c),  for  these  nonai spersive  solitons, 


s'  ~s  S’ 

Next,  we  use  the  basic  directional  form  [(7.8),  or  (8.5)],  now  for 
these  nondi spersive  solitons,  where  K<-  =  ^s/cs,  cf.  (7.17c).  The  result 
of  applying  (9.9)  therein  is,  finally  [3],  the  general  directional 
covariance 


Ks(Ar,t}  -  o^<  |B(a)|K1(|B[a)|  )> 


B(a)  =  “  [c  t  -  Ar  cos(A4>-a)], 

H  (9.10) 


with  Ar  =  r2_ri’  =  2 i »  as  before*  cf.  (8.5),  r  in  the  (x,y)-plane, 

of  course.  In  our  present  application  we  assume  semi -i sotropy ,  e.g., 

K<.  is  independent  of  direction  (in  the  downwind  regime),  so  that  we  have 


Atj)  =  a , 


(9.11) 


and  consequently  (9.10)  reduces  to  our  earlier  result  [3],  [4] 


9 

csx-Ar 

c  x-Ar 

\ 

Ks(Ar,T ) 

=  C 

semi -iso  °S 

°h/2 

Kl( 

oh/2 

(9.12) 
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This  reduces  at  once  to 


Ks(Ar,°)s. 


1  so 


2Ar 


2Ar 


K$(Ar ,0) , 


(9.13) 


(and  to  (9.7)  when  Ar  +  0). 

From  (9.13)  we  can  obtain  the  desired  wave  number  spectrum  V^2 (—  ’ *^ ) s ’ 
(2.3).  Noting  the  isotropy  of  our  model  (9.13),  we  have 


W„(k,0)  =  //  K  (Ar.Oe1-*—  d(Ar) 

2  S  [at]  S  ~  ~~ 

oo 

=  4Tta2a^  /  JQ(kAr)K^(2Ar/a^)(Ar)2  d(Ar), 
which  from  [30],  p.  410,  Eq.  (2)  becomes  finally 


(9.14a) 

(9.14b) 


w2(k,o)s  = 


2  2 
TT  0^0, 

S  h 


*C  l+(  ka.  /  2)3 


2  2  ^2  ^  ^  )  < 


(9.15) 


[In  the  present  application  k  =  2kQsin0oj,  (4oy  =  tt/2 ) ,  cf.  Sec.  2 
and  Fig.  2.1.] 

Finally,  we  comment  on  various  assumptions  underlying  our  present 
model.  We  note  that 

(1) .  The  soliton  speed,  c,,  is  treated  as  essentially  constant  here 

for  the  small  layer  thickness,  heff’  deduced  to  be  o(l-3  mm), 
cf.  Table  3.1,  this  means  that  h  ^  does  not  change  appreciably 
with  sizeable  changes  in  vQ:  changes  in  vQ  produce  changes  in  LQ, 
cf.  (9.3),  but  not  in  h^^.  This  latter  is  consistent  with  (9.5). 
Thus,  although  v^  can  have  sizeable  variances,  cf.  Sec.  (8.1) 
and  Eq.  (4.  2),  the  effective  layer  hg^  remains  unchanged,  once 
the  soliton  process  has  been  established.  While  this  appears 
plausible  and  gives  results  quite  consistent  with  the  experi¬ 
mental  data,  the  assumption  of  the  constancy  of  cs  needs  to  be 
directly  verified. 

(2) .  The  soliton  field  is  taken  to  be  (semi-)isotropic,  cf.  (9.11)  et 

seq.,  with  soliton  motion  in  the  direction  of  "downwind."  This 
semi-isotropicity  is  probably  not  often  the  actual  condition: 
the  soliton  field  is  anisotropic,  much  like  that  of  the  gravity- 
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capillary  wave  surface.  However,  we  expect  that  numerical  dif¬ 
ferences  are  subsumed  (i.e.,  disguised)  in  the  various  scaling 
parameters,  o^,  and  are  not  distinguishable  at  the  present 
level  of  experimental  refinement.  Again,  these  effects  need  to 
be  isolated  in  properly  controlled  experiments. 

The  detailed  structure  of  the  present  model  continues  to  be  investigated 
by  the  author. 

From  (9.9)  and  (9.15)  we  observe  the  limiting  behavior  of  these 
spectra  as  frequency  increases.  Thus,  we  have  (as  f  and  k -*-«>): 


W  (f  )  ~  8itc2 
s  s  s 


W2(k,0)s~  16tti 


(9.16) 


-4  -3 

showing  the  expected  k  and  f  structure.  Both  spectra  exhibit  this 
behavior  in  excellent  agreement  with  experiment  ([8],  [26]  for  k  - 
dependence,  for  example,  and  Fig.-  16  of  [32],  and  in  [13];  and  [14]  • 

_3 

for  f  -dependence).  [Of  course,  ultimately  these  dependences  must 
fall  off  faster,  in  order  to  ensure  the  physically  required  finite 

_A 

rms  slopes.  This  implies  a  fall-off  at  least  o(k  or  f  ),  or  if 

-4-f  ^  -3-f 

fractional  exponents  are  allowed,  at  least  as  0(k  or  f  ,  e>0).] 


With  the  results  of  Secs.  (8.1)-(8.3)  we  are  ready  to  determine 
the  spectral  spread,  AF,  whose  precise  form  we  define  below.  To 
evaluate  AF  we  return  to  (6.3)  and  insert  (8.3),  (8.19),  (8.29a)  as 
indicated.  The  result  is  explicitly. 


M(0)(t 

RT  lT 


(Sos2  ko9T9RT)NGS)(,ioi>e 


2ik(£ol*y0DT 


2  2„ 

-°at  n 


2  i  k 


o-ol 


■  Ar 


Ks(Ar,-r)d(Ar) , 


•  // 
[Ar] 


e 


(10.1) 
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where  now 


2  2  2  2 

o  =  0 oV  +  CO  R-  +  Rr 
A  o  6  o  ?  G 


(10.2) 


is  the  spread-coefficient.  Applying  (10.1)  to  (6.2)  gives  directly  the 
intensity  spectrum  of  the  received  (incoherent)  component  of  the  back- 
scatter,  X-<)0: 


W 


X-<X> 


2  2  — 
i  t-[\  00  -a.x  /2-icox  -iw  r+i2k  a_ ,  *v  nx 

(f)  =  2kM^B  /  dr  e  A  Re{Ko(T)ine  0  °~oi  oD 

-00 


2i k  a  ,  *Ar 


//  -ol- 

[Ar]  5 


d(Ar)  } 


(10.3) 


At  this  point  we  use  the  directional  form  (7.8)  for  the  soliton 
surface  waves,  e.g. , 

oo 

=  /  MO  <cos(co  x-K  -Ar  >  df  . 

3  o  0  b 

Inserting  (10.4)  into  the  integral  over  Ar  in  (10.3),  in  the  manner  of 
(7.9),  gives  us  for  these  solitons 


(10.4) 


\(VT)  =  ix 


i  k  I  2a  |  c  x 

'MW. e  5  • 


(10.5) 


which  reveals  the  Bragg  scatter  doppler  explicitly.  With  (10.5)  we  get 
for  (10.3)  at  once 


*W<f> = 


c  r 
0  1 


2  2  . 

-1C0X 


Re{K  (x) .  e 
o  'in 


-i  (co  +6co  )t 

0  0  }dx  , 


(10.6) 


where  the  mean  doppler  shift  is,  as  expected. 
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Sco 

o 


[2a  ,*v  n/c  +  2 [a  ,|c  /c  Ito 
L  — oi  -"OD'  o  u~oll  s  o  o 


(10.6a) 


cf.  (7.21),  and 


C  e  2k2G^^? 
o  o  2v 


Bt 


wJc  v„)w,  (4>  ) 

S'  s  o'  lVTvQ  a 


B  5  Rl  s2  koVR  2 


’  NGS^ojl’ 


(10.6b) 


with 


03  /s  /s 

v  =  2k  la  ,  |/2tt  =  -  {i  cos<f>  xsin9  T  +  i  sind)  Tsin0  T} ,  (10.6c) 

o  ou~oll  ttcq  -x  yoT  oT  ~y  yoT  oT 

cf.  (7.7).  (An  interesting  feature  to  note  about  our  result  (10.6) 
is  that  the  spectral  density  ^  of  the  soliton  surface  appears  as  a 
constant  factor  of  the  integral,  cf.  (10.8)  ff.) 

Observing  that  K  (x).jn  is  even,  cf.  (6.2b),  and  considering  only 
positive  frequencies  (f>0),  we  see  that  (10.6)  reduces  to 


(10.7) 


Equation  (10.7)  shows  the  spreading  and  frequency  shift  produced  in  the 
incident  (narrow-band)  signal  by  interaction  with  the  complex  wave 
surface  of  our  model.  If  long  cw  pulses  are  used,  i.e.,  long  compared 
to  the  decay  time  x.  =  a”*,  which  is  the  case  here,  we  can  replace 

KoCOinby  A^/2  for  an  effectively  cw  signal.  Then,  (10.7)  reduces  fur¬ 
ther  to 

C  A2  °°  -a2 t2/2 

^X-<X>^f)  =  2  e  cos(w-boc)x  dx  .  (10.8a) 

C  A2  r -  -b2  /2 a2 
_  oo  ,/n  oc'  A 

°A 


(10.8b) 
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The  (two-sided)  doppler  spectral  spread,  AF,  about  the  doppler 

shifted  carrier,  f  +Sf  ,  is  defined  by 
oo  J 


(10.9) 


where  the  angular  frequency  b  is  that  for  which  the  following  rela¬ 
tion  holds: 


^X-<X> 

^X-<X> 


(b  ) 
v  oc' 

lor- 


C  =  10  log1Qe  =  4.34. 


(10.10) 


Accordingly,  using  (10.8b)  in  (10.10)  gives  at  once 

bQC  »  log  C  =  1.714  aA;  (  \jz  log  t  =  1.714). 

Finally,  putting  (10.2)  together  with  (10.9)-(10.11)  we  obtain  the 
desired  expression  for  the  doppler  (frequency)  spread: 


(10.11) 


(10.12) 


Given  AF  by  experiment  at  two  different  frequencies  enables  us  to  esti- 
2 

mate  o0  and  R„,  as  we  describe  in  Sec.  4.  We  note  again  that  to  obtain  AF 

P  b 

here  it  is  not  necessary  to  evaluate  K^(At,t)  explicitly,  or  correspond¬ 
ingly,  cf.  (10.5).  The  same  is  true  of  the  doppler  shift  <SfQ,  cf. 

Sec.  7.  Of  course,  if  we  wish  the  spectral  density  W _ <  ,  or  the  back- 

scatter  cross  section  (cf.  Sec.  2),  the  structure  of  is  required. 
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